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Abstract 


Pulsed carbon dioxide (CO2) lasers have many applications in aeronautics, space research, 
weather monitoring and other areas. Full exploitation of the potential of these lasers is hampered 
by the dissociation of CO2 that occurs during laser operation. The development of closed-cycle 
CO2 lasers require active COO2 recombination (CO oxidation) catalyst and design methods for 
implementation of catalysts in CO2 laser systems. 

A monolith catalyst section model and associated design computer program, LAS CAT, are 
presented to assist in the design of a monolith catalyst section of a closed cycle CO2 laser system. 
Using LAS CAT, the designer is able to specify a number of system parameters and determine the 
monolith section performance. Trade-offs between the catalyst activity, catalyst dimensions, 
monolith dimensions, pressure drop, O2 conversion, and other variables can be explored and 
adjusted to meet system design specifications. 

An introduction describes a typical closed-cycle CO2 system, and indicates some 
advantages of a closed-cycle laser system over an open cycle system and some advantages of a 
monolith support over other types of supports. The development and use of a monolith catalyst 
model is presented. The results of a design study and a discussion of general design rules are 
given. Detailed instructions for the use of LAS CAT is given in Appendices 1 - 4 . 



Monolith Catalyst Model Development and Implementation 


L INTRODUCTION 

The design of a CO2 laser system is made under constraints of energy consumption, operating cost, 
size, weight, pressure drop, loop O2 conversion, catalyst activity, flow channel dimensions, particulate 
production, reliability, ease of operation, safety, and numerous others. Design decisions, made to satisfy 
one constraint, often push the limits imposed on other constraints. A means to design the monolith 
catalyst section of a closed-cycle CO2 laser system under selected restraint conditions is presented. 

A pulsed CO2 laser produces useful laser light output when the laser gas volume is exposed to a high 
electrical potential for a period on the order of milliseconds Low lying vibrational-rotational energy 
level CO2 molecules are excited to higher energy levels. The relaxation of the excited molecules to low 
lying energy levels produces infrared radiation (= 10 (im). The energy deposited by the electrical discharge 
alters the composition of the gas mixture by dissociation of CO2 into stoichiometric ratios of CO and O2. 
In a high pulse repetition rate (PRR) CO2 laser system, without a method of recombining the CO and O2, 
the laser gas O2 concentration would continue to increase with time. Depending on the particular system 
used, O2 concentrations above a threshold level, on the order of 1 %, severely degrades the laser output 
quality. The lasing process also raises the temperature of the laser gas volume. The output of the laser is 
sensitive to the temperature of the operating gas due to the effect of temperature on the distribution of CO2 
molecules among the rotational-vibrational energy levels. 

An open-cycle laser system requires a continuous supply of fresh gas to maintain O2 concentration 
below threshold levels. The higher the PRR, the higher the fresh gas flow rate required to maintain 
acceptable O2 concentrations. If rare isotope CO2 gases 1 ’ 3 , such as 12 C 18 02 , were to lx used in an open- 
cycle laser system, the cost of supplying large quantities of CO2 gases would be prohibitive. The safety 
requirements for handling and for the disposal of CO, O2, and CO2 gases is another concern for open- 
cycle systems. Clearly, for portable uses, such as in satellites or in the field, open-cycle system operation 
is not feasible. 

Closed-cycle systems can be envisioned to require an one time charge of the operating gas and to operate 
for a given number of pulses (> 10 7 ) at a specified PRR and power output 4 ' 7 . Portable closed-cycle 
systems would necessarily have power consumption constraints and, therefore, be limited in the output 
power and PRR. In a typical closed-cycle CO2 laser system (Fig. 1 ), a gas mixture, comprised of CO2 
and other "inert” gases, is continuously recirculated through the system by a blower. The recombination 
of CO and O2 is accomplished by the monolith catalyst section of the system (Fig. 2 ) and provides the 
laser section with a fresh supply of operating gas. The heat exchangers allow for the operation of the 
monolith catalyst section at elevated temperatures while maintaining moderate laser section temperatures. 
The monolith catalyst material is very porous (high surface area) and has catalytic material dispersed 
throughout. The performance of the monolith catalyst is measured by the section's ability to recombine 
CO and O2, the size and weight of the monolith, and the pressure drop produced as a result of gas flow 
through the monolith. Monolith catalyst performance is dependent on a number of interrelated factors, 
such as, the catalyst's geometry, convective heat and mass transport rates from the bulk gas to the catalyst, 
inlet gas properties, inlet molar gas flow rate, and catalytic activity. A monolith support is chosen over 
other supports, e.g. powders, beads, etc., because of the reduced level of particulate production, and the 
sturdiness of the monolith under high volumetric gas flow conditions. 
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Figure 1. A typical closed-cycle CO2 laser. 


2. MONOLITH CATALYST SECTION MODEL 

A flexible model of the monolith catalyst may be developed to determine the bulk-average gas 
temperature, composition, and pressure along the length of the monolith. The adjustable parameters 
required to specify the operating condition of the monolith are listed in Table 1. The model assumptions, 
balance equations, and model results follow. 



2.1 Model Assumptions 
•Steady state conditions. 

•Identical conditions exist in each channel of the monolith. 
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•Channel gas flow is laminar and fully developed. The neglect of entrance effects for heat and mass 
transfer calculations provide for a conservative estimate of the amount of catalyst required. Pressure drop 
due to entrance and exit effects is neglected, therefore the actual pressure drops will be somewhat greater 
than calculated. 

•The kinetics of the reaction, V 2 O 2 + CO CO 2 , is first order overall, reaction rate « [O 2 ] a [CO] b 
where sum of a and b equals one. O 2 and CO 2 appear in stoichiometric ratios, 2:1s C0:02. The 
concentration of O 2 is followed and the remaining species concentrations are calculated. The effect of 
temperature on the reaction rate is through an Arrhenius dependence of a reaction rate constant 
•Inlet gas composition, temperature and flow rate are known. Gas properties, viscosity, diffusivity, and 
thermal conductivity, are updated with changes in temperature, pressure, and reactant concentrations. 
•Diffusivity in the porous catalytic layer is calculated using size and void fraction of micropores and 
macropores. Equimolar counterdiffusion in porous catalytic layer is assumed. 

•Either adiabatic or isothermal monolith operating conditions can be selected by the operator. These two 
operating conditions provide the upper and lower bounds for oxygen conversion for a given inlet gas 
condition. 

• Axial heat conduction in the porous catalytic layer and support is assumed to be negligible. For adiabatic 
monolith operation, transverse heat conduction the porous catalytic layer and support is assumed to be 
such that the porous catalytic layer and support temperature is uniform transversely. 

•Slab geometry is used for calculation of the species concentration in the porous catalytic layer. A 
characteristic porous catalytic layer thickness is calculated to account for porous catalytic material in the 
channel comers. 

•Heat and mass transport between the flowing gas and the channel walls are described using the limiting 
Sherwood (Shoo) and Nusselt (Nuoo) numbers for constant wall concentration and temperature boundary 
conditions in square channels 8 . 


2.2 Model balance equations 

The steady state conservation equations for O 2 in the flowing gas and O 2 in the porous catalytic layer are 


dC 



, and 


dV 


dX 


2 

(p ¥ 


with initial and boundary conditions, 

O g = 1 at C = 0 , ¥=1 at X = 0 for all £ , and ^ 

dX 


( 1 , 2 ) 

0 at X = 1 . (3,4,5) 


Og is the dimensionless bulk- average O 2 molar flow rate in the flowing gas. Variables and parameters are 
defined in detail in the notation list. Ow is the O 2 concentration at the channel wall times bulk-average 
volumetric flow rate divided by the inlet bulk-average O 2 molar flow rate. [ O g -<t>w ] is proportional to 
the concentration driving force for transport of O 2 from the flowing gas to the channel wall. £ is the 
dimensionless distance down the length of the monolith. ¥ is the dimensionless O 2 concentration inside 
the porous catalytic layer at a dimensionless depth X into the layer. The solution to equation (2) yields an 
overall O 2 reaction rate in the porous catalytic layer as a function of Ow- At steady state, this reaction rate 
is equal to the rate of transport of O 2 from the flowing gas to the channel wall, and leads to 


Ow = C>g 


_L_ 
.oc + 1. 


( 6 ) 


The steady state energy balance on the flowing gas relates the rise in the dimensionless bulk-average 
temperature of the flowing gas, 0g, to the heat transferred to the gas from the channel wail and yields 
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, with initial condition, 0 G = 1 at £ = 0 , (7,8) 
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where 9 g is the dimensionless bulk-average temperature of the flowing gas, and 9w is the dimensionless 
temperature at the channel wall. An energy balance on the porous catalytic layer equates the heat 
generation from the oxygen consumption reaction to the heat transferred from the porous catalytic layer 
and yields 


OLfiL 


dty = [Og - 9w] 


(9) 


Combining the two energy balance equations, (7) and (9), and substituting for O w using equation (6) 
results in two final differential equations, to be integrated. 


dO G _ 

dC 
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an d = [dhSgmS ; ] ®o 
[ 1 + a] P 


( 10 , 11 ) 


with the same initial conditions for O g and 0 G as above. Pressure drop for the laminar flow is calculated 
using the Hagen-Poiseuille equation. 


— — = „ , with an initial condition of p = l at £ = 0 , (12,13) 

d£ Eu Re 

where p is the dimensionless pressure. Equations (10), (11), and (12) are coupled to each other through 
the variable O g and the parameters a, to, T, y, St, Eu, and Re, each of which is a function of 0 G , <X> G and 

P. 

2.3 Model results 

The adjustable parameters required to specify the monolith catalyst section operating conditions are listed 
below in Table 1. They are used to compute the dimensionless parameters a, co, T, y, St, Eu, and Re, and 
to integrate the dimensionless variables in equations (10), (11), and (12). A computer program (see 
section 5 ) was written to perform integration using fourth order Runge-Kutta method. The program 
outputs 0 G , 9w, O g , Ow and p as C varies from zero to the desired dimensionless monolith catalyst 
section length. The computer program requires readily available parameters to calculate parameters such 
as mass and heat transfer coefficients, bulk-gas and effective diffusion coefficients, and thermal 
conductivity. Results for a monolith catalyst section operating under conditions specified in Table 1 are 
shown in Figs. 3, 4 and 5. 

Fig. 3 shows the behavior of 0, 0 , and p as £ varies from 0 to 100. 0 G =0w and <f> G =Oy/, therefore 
0 G , 0w» ^g. and <I>w are not shown individually on this small scale plot. The behavior of 0, O, and p 
can be explained using equations (1)-(12), in which all parameters and variables are positive valued. The 
right side of equation (11) is positive, therefore, d0 G /d£ > 0. The left hand side of equation (9) is 
negative, therefore, 0w > 0 g- Similar examination of equations (10) and (6) results in dOc/d^ < 0 and 
O w > ^g- Equation (12) indicates that dp/dC, < 0. The drop in p is negligible for the gas flow rate 
specified in Table 1. Pressure drop can be significant for higher gas flow rates. 

Fig. 4 is an expansion of Fig. 3 in the £ = 45 to 50 region. The separation between 0 G and 0w is now 
visible and provides an indication of the thermal driving force between the channel wall and the bulk-gas. 
Heat produced from the oxidation of CO in the active catalyst layer of the channel wall is transferred from 
the channel wall to the bulk-gas. 

Fig. 5 is also an expansion of Fig. 3 in the £ = 45 to 50 region. An additional curve, v F(A.=l)*Ow, is 
the scaled O 2 concentration at the centerline of the monolith channel support wall. O 2 is consumed in the 
porous catalytic layer as mass transport to the center of the monolith channel support wall occurs, 
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Figure 6. Optimum design vs. standard commercial size. 







therefore v F(X=l)*<I>w is less than <X>w or <I>g. The separation between the curves indicates the O 2 mass 
transport driving force. The difference between C>g and d>w is an indication of the driving force between 
the bulk gas and the channel wall. The difference between Owand4 / (X?=l)*<I>wis an indication of the 
driving force between the channel wall and the centerline of the monolith channel support wall. 


Table 1 : Monolith catalyst section operating parameters 


Monolith dimension 
Facial cross sectional area = 1.0e+4 mm 2 
Support wall thicknessfincluding active layer) = 1mm 
Active layer thickness = 0.25mm 
Gas inlet properties 
Flowrate = 0.25 liters/sec 
Temperature = 300 K 
Pressure = 101.325 kPa 
Composition 37% CO 2 , 2% CO, 1% O 2 
40% N 2 , 20% He 
Thermal operation is adiabatic. 


Catalyst properties 
Void fraction as macropores = 0.25 

Void fraction as micropores = 0.48 
Macropore radius = 500nm 
Micropore radius = 12nm 
Active layer density = 0.5e-3 g/mm 3 
Reaction rate constant = 123.4mm 3 /gcat sec 
(at 298K) 

Activation energy constant = 39700 J/mol 


3. DESIGN STUDY 


3 J Des i gn co n s tr aints 

A design study for a monolith catalyst section, operating under parameters similar to those listed in 
Table 1, was performed. For the study, the support wall was assumed to be composed entirely of active 
catalytic material. Additional constraints, 25% O 2 conversion (Oc[exit]=0.75) and a 0.125kPa (=0.5 in 
H 2 O) pressure drop across the monolith section, were imposed. The monolith facial cross sectional area 
and inlet gas conditions (composition, temperature, flow, and pressure) were held constant. Active 
support wall thickness and channel opening dimensions were varied and the minimum monolith length 
determined under the imposed constraint conditions. 

The rationale for minimizing the monolith length follows from the assertion that the smallest monolith 
leads to the smallest, lightest and least expensive laser system. In portable laser systems size, weight, and 
cost considerations are critical. The choice of the shortest monolith section reduces the system weight by 
reducing the required monolith section housing length. The housing material weight per unit length is 
typically an order of magnitude more than that of the monolith material. The monolith material can be 
quite expensive; a smaller system requires less materials and is less costly. The 25% conversion 
requirement ensures a constant gas temperature rise. A constant pressure drop across the monolith section 
ensures that an identical amount of energy is expended to circulate the gas through the system. Excluding 
laser pulse energy, system energy requirements (blower energy and heating and cooling loads) are fixed 
by the pressure drop and conversion requirements. Constant gas inlet conditions and 25% conversion 
ensure that a chosen laser PRR can be maintained without exceeding a the maximum allowable laser 
section inlet O 2 concentration. 

3.2 Design results 

Fig. 6 shows the optimum monolith section design geometry in comparison with a standard monolith 
design. The optimum monolith length is substantially less than the standard monolith length, whereas the 
optimum support wall thickness is substantially greater than standard thickness. As the monolith support 
wall thickness is increased, more catalytic material can be packed into a shorter monolith while conversion 
and pressure drop constraints are still satisfied. However, a thicker support wall has a larger O 2 mass 
transport resistance. A point is reached where the benefit of having thicker walls is negated by the large 
O 2 mass transport resistance and an optimum monolith length is determined. 
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4, DI SCUSSION 


The design study presented above shows that the use of off-the-shelf commercial monolith designs for 
CO 2 laser applications dramatically increases the overall monolith size required relative to the optimum 
design presented here. Off-the-shelf monoliths are designed for pollutant emission control and have been 
optimized for large gas flows and fast reactions at high temperatures. They have relatively thin monolith 
catalyst section dimensions and high % void volume. The slow reaction rates obtained over laser catalysts 
allow use of relatively thick monolith wall dimensions. With lower gas flow rates, low % void volume 
can be used to obtain compact monoliths, thereby reducing the size, weight, and cost of the laser system. 

Using a computer program to generate variable and parameter values along the length of the monolith 
allows for rapid optimization of monolith section under a set of constraints. Use of a computer program 
also allows for complex channel geometries, such as cylindrical, hexagonal, triangular, etc., to be 
incorporated into the monolith design. 


5. MODEL IMPLEMENTATION USING LASCAT 

LASCAT, a computer program based on the monolith catalyst section model presented above, provides 
a means to design a monolith catalyst section that will satisfy a user specified set of design requirements. 
LASCAT requires the specification of the parameters listed in Table 1 and an few others. Values of key 
parameters ( Re, Cpg, pg, Dab. DABeff. and mole fractions ) are provided at the inlet and exit of the 
monolith section. Fq, T and P are integrated instead of the corresponding nondimen sional parameters, 
$G, 6g and p,C,T and P values are provided along the length monolith. See Appendix 2, a step by step 
tutorial, for an example of the LASCAT output . 

Appendix 1 gives detailed information on the implementation and compatibility of LASCAT. Appendix 3 
details each menu operation of LASCAT. Appendix 4 is a program listing of LASCAT and provides a 
wealth of information on program structure, variable definitions and units, and physical properties 
computational methods and references. 


6. NOTATION LIST 


6.1 Variables and parameters 
A = CO oxidation reaction rate constant ( mm 3 / g-cat s) 

C = concentration of oxygen (mol / mm 3 ) 

Cpc = bulk-average gas heat capacity ( J / g K) 

Dab = diffusion constant for O 2 in bulk gas mixture ( mm 2 / s ) 

^ABeff = effective diffusion constant for oxygen in 
porous catalytic layer ( mm 2 / s ) 
dh = hydraulic diameter of monolith channel, 4(cross sectional area )/( wetted perimeter) ( mm ) 
Eu = Euler number, Eu = Po dh 4 / d 2 p g 

F g = bulk-average gas molar flow rate of oxygen, F G = Cgd ( mol / s ) 
h = average heat transfer coefficient ( W / mm 2 K ) 

AH^rz heat of CO oxidation reaction ( J / mol O 2 convened) 
k = bulk-average gas thermal conductivity (J / m K ) 
km = average mass transfer coefficient ( mm / s ) 

Nuoo = limiting Nusselt number, Nuoo = h dh / k 
P = pressure in the channel at a position C, (kPa) 
p = dimensionless pressure at position £, p = P / Po 
Re = Reynolds number , Re = dh D p g / (igdh 2 

S = Surface area of the channel wall per unit open volume of monolith channel ( mm 2 / mm 3 ) 
Shoc= limiting Sherwood number, Shoo= k m dh / Dab 
S t = Stanton number , St = h dh 2 / Do p g Cp G 
tc = characteristic thickness of porous catalytic layer ( mm ) 


6.2 Subscripts 

C,c in porous catalytic layer 

G,g in flowing gas 

0,o at inlet conditions 

W,w in gas at channel wall 
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T = temperature ( K ) 

x = distance along the length of the monolith channel ( mm ) 
z = depth into the porous catalytic layer ( mm ) 


Greek 

a = dimensionless rate constant, a = tcP c T|A / k m 

y= dimensionless mass transfer coefficient, y= dh S k m dh 2 / t>o 

r = dimensionless volumetric flow rate, T = d / Do 

co = dimensionless heat of reaction, CD = -AH rxn k m FG 0 / h T G0 Do 

C, = dimensionless distance along the length of the monolith channel, £ = x / dh 

tj = effectiveness factor of the reaction in the monolith porous catalytic layer , t\ = tanh (cp) / cp 

6 = dimensionless bulk gas temperature, 0 g=Tg / Too, 0w=Tw/ Too 

X - dimensionless depth into the porous catalytic layer, X = z / tc 

|lg = viscosity of gas in channel ( Pa-s ) 

pc = density of porous catalytic layer material ( g-cat / mm 3 ) 

pg = density of flowing gas ( g / mm 3 ) 

d = volumetric flow rate of flowing gas ( mm 3 / s) 

d>G = dimensionless bulk-average oxygen concentration in flowing gas, O g = F G / Fgo 
d>w= dimensionless channel wall surface oxygen concentration, <t>w= [Cw(0 d(£)] / Fgo 
< p = Thiele modulus, cp 2 = [t c 2 p c ‘n A 3 / DABeff 

'F = dimensionless O 2 concentration in the porous catalytic layer. 'F(X) = [Cc(^) /C W (C)] 
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Appendix 1: Program Implementation and Compatibility 


The LASCAT program is written in the FORTRAN programming language and is compatible 
with FORTRAN 77 standards. Detailed implementation instructions are provided below for 
Apple Macintosh and Digital Equipment Corporation VAX computers. Familiarity with the 
particular operating system is assumed. 


Computer: Apple Macintosh SE or II. 

Operating System: Apple System Version 6.0 

Application Program: Absoft's MacFortran/020 Version 2.3 

Compilation Options (Macintosh SE): 

-Use compilation options B (Compile Using Long Addresses) 

E (Generate Errors List), and 
U ( * = Unit 9 ). 

Compilation Options (Macintosh II): 

-Use compilation options B (Compile Using Long Addresses), 

E (Generate Errors List), 

M (68020/68030 instructions), 

P (68881/68882 instructions), and 
U (* = Unit 9). 

Compilation and Execution: 

-Place the LASCAT program in the same folder as the MacFortran/020 application and 
supporting files. 

-Double click on MacFortran/020 application to launch the MacFortran/020 application. 

-From the File Menu choose the "Select File" option. Select File LASCAT. 

-From the Compile Menu choose the "Options" option and verify the compilation options 
selected are identical to the compilation options described above. If the options selected 
need to be changed, remember to save the new set of compilation options by clicking on the 
"Save" box. Exit the compilation options section by clicking on the "OK" box. 

-From the Compile Menu choose the "Compile and Execute" option. The program will be 
compiled and executed. The output file, LDATA, and any selectable parameter files will be 
placed in the same folder as the LASCAT and MacFortran/020 files. 


Computer: Digital Equipment Corpo ration VAX 11/780 
Operating System: VAX/VMS Version 4.6 
Application Program: VAX FORTRAN V4.8-276 
Compilation and Execution : 

-Rename the LASCAT program LASCAT.for 
-To compile LASCAT.for, type for LASCAT.for 
-To link LASCAT.for type link LASCAT 
-To execute LASCAT.for type run LASCAT 

-The output file, LDATA.DAT, and any selectable parameter files will be placed in the your 
directory. Use an editing program to examine files. 
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Appendix 2: LASCAT Tutorial 


The LASCAT Tutorial provides a step by step example of the use of the program LASCAT. 
Required user input and comments are detailed on the left margin and the corresponding screen 
outputs are indented. Input(s) required by the operater are boldfaced. The example shown 
below was run on a Macintosh II using Absoft's MacFortran/020 Version 2.3 (Review the 
Program Implementation and Compatibility section for program complilation and execution 
details.). This tutorial will be of the greatest benefit if used while actually executing the 
program. 

Compile and execute LASCAT application program. 

...... PROGRAM LASCAT ******* 

The purpose of this program is to calculate the gas concentration and 
temperature profiles of a monolith catalyst section of a C02 laser. The 
C02 decomposes when the laser is pulsed. The CO and 02 produced as a result 
of pulsing are detrimental to the efficient operation of the laser. The 
recombination reaction is CO + 1/2 02 -> CO. This program provides the 
means to model the performance of a monolith catalyst section under various 
gas compositions, temperatures, catalyst activities, gas flowrates, oxygen 
conversion, monolith face and length dimensions. Results can indicate if 
constraints such as conversion, maximum gas temperature, monolith weight 
are satisfied and how the system parameters may be altered to meet these 
constraints. Parameters and options may be altered to taylor the monolith 
design. Default values can also be used as a starting point for the design 
process. A review of the parameters and options chosen may be made prior 
to execution of the computational portion of the program. 

(HIT RETURN TO CONTINUE) 

A program introductory statement is presented. Type RETURN (key) after reading. 

LASCAT Main Menu 

1 ) Read in new operating parameters 

2) Show current operating parameters 

3) Change operating parameters 

4) Run program 

5) Exit program 

Type in number corresponding to choice above. 

The LASCAT Main Menu is presented. The main menu provides different options (1-5). 
Specific details for each option are provided in the Menu Description section. To choose one of 
these options, simply type in the desired number, then RETURN. For this example, type 2 
and RETURN. 


SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 
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Channel inner dimension: 4.00 No. Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flow rate(liters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 

Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): 02 conversion. %: 2.500 

Computation loop step size(mm): .5556 Display every 5.00 mm 

==> Hit Return when finished viewing <=*= 

The selectable parameter summary lists the parameter values that will be used in the 
computation portion of the program. When the program is initially run, default parameter 
values are assigned. As we'll see later, the values of the parameters can be changed. Hopefully, 
a review of the parameter summary above will hint to the meaning of each parameter. The use 
of each parameter is detailed in the Menu Description section. Let's assume the present values 
are satisfactory and proceed. Type RETURN. 

LASCAT Main Menu 

1 ) Read in new operating parameters 

2) Show current operating parameters 

3) Change operating parameters 

4) Run program 

5) Exit program 

Type in number corresponding to choice above. 

We have returned to the main menu. To save space, the main menu will listing will be 
abbreviated as "LASCAT Main Menu....". To run the computational portion of the program using 
the current parameter set listed in the selectable parameter summary type 4 and RETURN. 
Note: the program will alert vou that it has completed computations bv beeping three times . 

*****SEE FILE LDATA FOR RESULTS* — 


Initial values 

Reynolds number = 9.880 

Gas Heat Capacity (J/K-g) = 1 .797 

Gas Density (g/cm A 3) « .001 1 81 7 

Gas Velocity (mm/s) = 39.063 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 
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Effective Diffusivity(cm A 2/s)= .053407 

Step size(mm) = .555556 


Distance 

%Conver 

02gas 02wall 

02center 

Tgas 

Twall 

DPress 

(mm) 


<SBSS 

=(mMOLA)===> 

(Kelvin) 

(Kelvin) 

(kPa) 

0.000 

0.0000 

.4065 

.4060 

.4058 

300.000 

300.179 

-0.000E+00 

5.000 

0.2239 

.4052 

.4048 

.4046 

300.236 

300.416 

-0.730E-05 

10.000 

0.4574 

.4040 

.4035 

.4033 

300.474 

300.656 

-0.1 46E-04 

15.000 

0.6933 

.4027 

.4022 

.4020 

300.714 

300.897 

-0.2 1 9E-04 

20.000 

0.9314 

.4014 

.4009 

.4008 

300.956 

301.141 

-0.293E-04 

25.000 

1.1718 

.4001 

.3996 

.3994 

301.201 

301.388 

-0.366E-04 

30.000 

1.4146 

.3988 

.3983 

.3981 

301.448 

301.637 

-0.440E-04 

35.000 

1.6598 

.3975 

.3970 

.3968 

301.697 

301.888 

-0.51 3E-04 

40.000 

1 .9075 

.3961 

.3957 

.3955 

301.949 

302.141 

-0.587E-04 

45.000 

2.1577 

.3948 

.3943 

.3941 

302.204 

302.398 

-0.661E-04 

50.000 

2.4104 

.3934 

.3930 

.3928 

302.461 

302.656 

-0.735E-04 

52.778 

2.5519 

.3927 

.3922 

.3920 

302.605 

302.801 

-0.776E-04 


Final values 
Mole fractions: 

He = .2001 Ar = 0.0000 C02 * .3706 
CO = .0195 02 = .0098 N2 = .4001 


Gas Pressure (kPa): 101.325 

Reynolds number = 9.817 

Gas Heat Capacity (J/K-g) = 1 .799 

Gas Density (g/cm A 3) = .0011720 

Gas Velocity (mm/s) = 39.378 

Effectivness factor = .99965 

Bulk Gas Diffusivity(cm A 2/s)= .22832 

Effective Diffusivity(cm A 2/s)= .053988 

Step size(mm) = .555556 


==> Hit Return to Return to Main Menu <== 

The information above provides initial and final values of important parameters along with 
values of key parameters along the length of the monolith. The units for each parameter are 
specified. A few minutes spent in reviewing the trend of each parameter, either initial vs. final 
or along the length of the monolith, is well worth the time. For example, the gas heat capacity 
increases from 1.797 to 1.799 due to the change in gas composition and the change in gas 
temperature. The gas temperature rises from 300.000 to 302.605 due to the exothermic 
nature of the reaction CO + I/2O2 -> CO2 . Type RETURN. 

LASCAT Main Menu.... 

To change the value of a parameter or parameters, type 3 and RETURN. 

Parameter Modification Menu 

1) Change monolith physical dimensions 

2) Change inlet gas composition 

3) Change inlet gas volumetric flow rate 

4) Change inlet gas temperature 

5) Change inlet gas pressure 

6) Change catalyst activation energy 
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7) Change catalyst reaction rate constant 

8) Change catalyst active layer density 

9) Change macro/micro pore radius & void-fraction 

10) Change thermal operation(Adiabatic/lsothermal) 

11) Change output profile (Full/ Summary) 

12) Change termination variable(Conversion/Length) 

13) Change computation display interval 

14) Show current operating parameters 

15) Save current operating parameters 

16) Run program 

17) Return to main menu 

18) Exit program 

Type in number corresponding to choice above 

You have entered the Parameter Modification Menu. The Parameter Modification Menu gives a 
number of options (1*18). The parameters are detailed in the Program Operations section. 
The Parameter Modification Menu will be abbreviated "Parameter Modification Menu....". Let's 
change the termination variable (option 12). Type 12 and RETURN. 

Integration can be stopped by either specifying 
an amount of 02 conversion (0), or by specifying 
a monolith length (L). 

Currently , integration will be stopped 
by specifying an amount of 02 conversion (0). 

Type the letter 0 or the letter L. Hit return. 

The letter L is for program computation section termination on monolith length. The letter 0 is 
for termination on O 2 conversion and requires a conversion length specification. Type L and 
RETURN. 


Current desired monolith length «= 5.000 (mm) 
input desired monolith length (5 - 4000mm) 

The current restriction on the desired monolith termination length is given. The restrictions 
presented must be followed. Attempts have been made to prevent out of bounds values from 
being input, however "nothing is perfect". Try entering 2 (this shorter than the minimum 
5mm ). The program will force a proper input. Type 100 and RETURN. 

Parameter Modification Menu.... 

Let’s see what has been changed. Type 14 and RETURN. 

SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No.Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 
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Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flow rate(liters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 

Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 
Change ==> Termination on (02 conversion/Length): Length. Length(mm): 100.000 
Computation loop step size(mm): .5556 Display every 5.00 mm 

==> Hit Return when finished viewing <== 

Type RETURN. To avoid changing from the default parameter values to frequently used values, 
a parameter set can be saved to a file. These files can be read from at a later time. Type 15 
and RETURN. 

Would you like to save the current operating 

parameters to a New(N) or Existing(E) file? 

Were creating a new file. Type N and RETURN. 

What's the name of the NEW 
parameter file to be opened? 

Name this file CATONE (abbreviation for CATALYST-ONE ). Type CATONE and RETURN. 

Parameter Modification Menu.... 

Using the revised parameter set, run the computational portion of the program by typing 16 
and RETURN. 


***’*SEE FILE LDATA FOR RESULTS***’ 
Initial values 

Reynolds number = 9.880 

Gas Heat Capacity (J/K-g) = 1 .797 

Gas Density (g/cm A 3) = .0011817 

Gas Velocity (mm/s) = 39.063 

Effectless factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 
Effective Diffusivity(cm A 2/s)= .053407 
Step size(mm) = .555556 


Distance 

(mm) 

0.000 


%Conver 02gas 02wall 02center Tgas 
<==s=(mMOL/L)==*=> (Kelvin) 
0.0000 .4065 .4060 .4058 300.000 


Twall DPress 

(Kelvin) (kPa) 

300.179 -0.000E+00 
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5.000 

0.2239 

.4052 

.4048 

.4046 

300.236 

300.416 

-0.730E-05 

10.000 

0.4574 

.4040 

.4035 

.4033 

300.474 

300.656 

-0 . 1 46E-04 

15.000 

0.6933 

.4027 

.4022 

.4020 

300.714 

300.897 

-0.21 9E-04 

20.000 

0.9314 

.4014 

.4009 

.4008 

300.956 

301.141 

-0.293E-04 

25.000 

1.1718 

.4001 

.3996 

.3994 

301.201 

301.388 

-0.366E-04 

30.000 

1.4146 

.3988 

.3983 

.3981 

301.448 

301.637 

-0.440E-04 

35.000 

1.6598 

.3975 

.3970 

.3968 

301.697 

301.888 

-0.51 3E-04 

40.000 

1.9075 

.3961 

.3957 

.3955 

301.949 

302.141 

-0.587E-04 

45.000 

2.1577 

.3948 

.3943 

.3941 

302.204 

302.398 

-0.661 E-04 

50.000 

2.4104 

.3934 

.3930 

.3928 

302.461 

302.656 

-0.735E-04 

55.000 

2.6657 

.3921 

.3916 

.3914 

302.720 

302.918 

-0.809E-04 

60.000 

2.9236 

.3907 

.3902 

.3900 

302.983 

303.182 

-0.883E-04 

65.000 

3.1842 

.3893 

.3888 

.3886 

303.247 

303.449 

-0.957E-04 

70.000 

3.4475 

.3879 

.3874 

.3872 

303.515 

303.718 

-0.1 03E-03 

75.000 

3.7136 

.3865 

.3860 

.3858 

303.785 

303.990 

-0.111 E-03 

80.000 

3.9825 

.3851 

.3846 

.3843 

304.058 

304.265 

-0.1 1 8E-03 

85.000 

4.2543 

.3836 

.3831 

.3829 

304.334 

304.543 

-0.1 26E-03 

90.000 

4.5290 

.3822 

.3817 

.3814 

304.613 

304.824 

-0.1 33E-03 

95.000 

4.8066 

.3807 

.3802 

.3800 

304.895 

305.108 

-0.141 E-03 

100.000 

5.0873 

.3792 

.3787 

.3785 

305.180 

305.395 

-0.1 48E-03 


Final values 
Mole fractions: 

He = .2001 Ar = 0.0000 C02 - .3712 
CO= .0190 02 - .0095 N2 - .4002 
Gas Pressure (kPa): 101.325 

Reynolds number - 9.753 

Gas Heat Capacity (J/K-g) = 1 .800 

Gas Density (g/cm A 3) = .0011624 

Gas Velocity (mm/s) = 39.693 

Effectless factor = .99961 

Bulk Gas Diffusivity(cm A 2/s)= .23125 
Effective Diffusivity(cm A 2/s)= .054566 
Step size(mm) = .555556 


— > Hit Return when finished viewing <== 

In addition to the screen output the selectable parameter summary and the section above are 
written to file LDATA(Macintosh) or LDATA.DAT(VAX). (For Macintosh only: Anv additional 
"runs" in this session will be tacked on to the end of LDATA. If the program is terminated and 
rerun the old LDATA file will be erased and a new LDATA file started. If vou desire to keep the 
information in LDATA. rename. LDATA prior to running LASCAT again . ) Type RETURN. 


LASCAT Main Menu.... 

After each computational run, you are returned to the main menu. Let's change another 
parameter. Type 3 and RETURN. 

Parameter Modification Menu.... 

To change the inlet gas volumetric flow rate, type 3 and RETURN. 


Current inlet volumetric 
flow rate - .250 (I/s) 
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Input monolith inlet volumetric flow rate 
The desired flowrate is 2.5 liters/sec. Type 2.5 and RETURN. 

Parameter Modification Menu.... 

Let’s examine the change to the parameter set. Type 14 and RETURN. 

SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No.Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 

Change -=> Gas Flow rate(liters/s): 2.500 Gas Temperature(K): 300.00 
Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 

Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): Length. Length(mm): 100.000 

Computation loop step size(mm): 5.0000 Display every 5.00 mm 

==> Hit Return when finished viewing <== 


Type RETURN. 


Parameter Modification Menu.... 

Run computational portion of the program by typing 16 and RETURN 

—••SEE FILE LDATA FOR RESULTS”*** 

Initial values 

Reynolds number = 98.802 

Gas Heat Capacity (J/K-g) = 1 .797 

Gas Density (g/cm A 3) = .0011817 

Gas Velocity (mm/s) = 390.625 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 

Effective Diffusivity(cm A 2/s)= .053407 

Step size(mm) = 5.000000 

Distance %Conver 02gas 02wall 02center Tgas Twall DPress 
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(mm) 


<== 

=(mMOL/L)===> 

(Kelvin) 

(Kelvin) 

(kPa) 

.000 

.0000 

.4065 

.4060 

.4058 

300.000 

300.179 

-0.000E + 00 

5.000 

.0156 

.4064 

.4059 

.4057 

300.024 

300.203 

-0.730E-04 

10.000 

.0391 

.4062 

.4058 

.4056 

300.047 

300.227 

-0.1 46E-03 

15.000 

.0624 

.4061 

.4057 

.4055 

300.071 

300.250 

-0.21 9E-03 

20.000 

.0855 

.4060 

.4055 

.4054 

300.095 

300.274 

-0.292E-03 

25.000 

.1086 

.4059 

.4054 

.4052 

300.119 

300.298 

-0.365E-03 

30.000 

.1317 

.4057 

.4053 

.4051 

300.142 

300.322 

-0.438E-03 

35.000 

.1549 

.4056 

.4052 

.4050 

300.166 

300.346 

-0.51 1 E-03 

40.000 

.1781 

.4055 

.4050 

.4049 

300.190 

300.370 

-0.584E-03 

45.000 

.2013 

.4054 

.4049 

.4047 

300.214 

300.394 

-0.657E-03 

50.000 

.2245 

.4052 

.4048 

.4046 

300.238 

300.418 

-0.730E-03 

55.000 

.2478 

.4051 

.4047 

.4045 

300.261 

300.442 

-0.803E-03 

60.000 

.2711 

.4050 

.4045 

.4043 

300.285 

300.466 

-0.877E-03 

65.000 

.2944 

.4049 

.4044 

.4042 

300.309 

300.490 

-0.950E-03 

70.000 

.3177 

.4047 

.4043 

.4041 

300.333 

300.514 

-0.1 02E-02 

75.000 

.3411 

.4046 

.4041 

.4040 

300.357 

300.538 

-0.1 10E-02 

80.000 

.3644 

.4045 

.4040 

.4038 

300.381 

300.562 

-0.1 17E-02 

85.000 

.3879 

.4043 

.4039 

.4037 

300.405 

300.587 

-0.1 24E-02 

90.000 

.4113 

.4042 

.4038 

.4036 

300.429 

300.611 

-0.1 32E-02 

95.000 

.4347 

.4041 

.4036 

.4034 

300.453 

300.635 

-0.1 39E-02 

100.000 

.4582 

.4040 

.4035 

.4033 

300.477 

300.659 

-0.1 46E-02 


==> Hit Return to Return to Main Menu <== 


Notice that the program computation terminated when the monolith length reached 100 mm, 
whereas, previously, the termination condition was an oxygen conversion of 2.5 %. Type 
RETURN. 


LASCAT Main Menu.... 

The parameter set has been changed from the CATONE set. Let's return to the CATONE parameter 
set by reading in the CATONE set file. Type 1 and RETURN. 

Read in data from existing parameter file?(Y/N) 

Type Y and RETURN. 

What is the name of the existing parameter file to be opened? 

Type CATONE and RETURN (Note that the required file name to be inputted may be lower 
/upper case sensitive). 


LASCAT Main Menu.... 

Let's review the parameter set. Type 2 and RETURN. 

SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No. Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 
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Change ==> 


Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flow rate (I iters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 

Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): Length. Length(mm): 100.000 

Computation loop step size(mm): .5556 Display every 5.00 mm 

==> Hit Return when finished viewing <== 


Type RETURN. 


LASCAT Main Menu.... 

Let's explore the use of some other parameters. Type 3 and RETURN. 

Parameter Modification Menu.... 


Type 11 and RETURN. 


Full concentration and temperature profile (P) or summary (S)? 
Type the letter P or the letter S. HIT RETURN. 

Choose summary by typing S and RETURN. 

Parameter Modification Menu.... 


Type 16 and RETURN. 

— *SEE FILE LDATA FOR RESULTS***** 


Initial values 

Reynolds number * 9.880 

Gas Heat Capacity (J/K-g) - 1.797 

Gas Density (g/cm A 3) « .0011817 

Gas Velocity (mm/s) = 39.063 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 

Effective Diffusivity(cm A 2/s)= .053407 

Step size(mm) = .555600 


Twall DPress 

(Kelvin) (kPa) 
300.179 -0.000E + 00 
305.395 -0.1 48E-03 


Distance %Conver 02gas 02wall 02center Tgas 
(mm) <===(mMOLA)===> (Kelvin) 

.000 0.0000 .4065 .4060 .4058 300.000 

100.008 5.0878 .3792 .3787 .3785 305.180 
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Final values 
Mole fractions: 

He = .2001 Ar * 0.0000 C02 = .3712 
CO = .0190 02 = .0095 N2 = .4002 


Gas Pressure (kPa): 101.325 

Reynolds number = 9.753 

Gas Heat Capacity (J/K-g) = 1.800 

Gas Density (g/cm A 3) = .0011624 

Gas Velocity (mm/s) = 39.693 

Effectivness factor = .99961 

Bulk Gas Diffusivity(cm A 2/s)= .23125 

Effective Diffusivity(cm A 2/s)= .054566 

Step size(mm) = .555600 


==> Hit Return to Return to Main Menu <== 

The key parameters along the length of the monolith are presented for the initial and final 
(termination condition met) monolth distances. The file LDATA also has this summary form. 

Type RETURN. 

LASCAT Main Menu.... 

Please feel free to play around with the program. To exit the program by type 5 and 
RETURN. Take a look at files, LDATA and CATONE. A listing of both files is given 
below. The table entries are separated by tab characters for further processing by programs 
such as Cricket Graph and Excel (these are both Macintosh programs). Adjusting the tabs to 1 
space may help to align the table. 

File LDATA 


SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No.Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flowrate(liters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 
Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): 02 conversion. %: 2.500 

Computation loop step size(mm): .5556 Display every 5.00 mm 


Initial values 

Reynolds number = 9.880 

Gas Heat Capacity (J/K-g) = 1.797 

Gas Density (g/cm A 3) = .0011817 

Gas Velocity (mm/s) = 39.063 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 
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Effective Diffusivity(cm A 2/s)= .053407 
Step size(mm) = .555556 


Distance 

%Conver 02gas 

02wall 02center 

Tgas 

Twall D Press 

(mm) 


<=== 

:(mMOL/L)=i 

==> 

(Kelvin) 

(Kelvin) (kPa) 

0.000 

0.0000 

.4065 

.4060 

.4058 

300.000 

300.179 -0.000E+00 

5.000 

0.2239 

.4052 

.4048 

.4046 

300.236 

300.416 -0.730E-05 

10.000 

0.4574 

.4040 

.4035 

.4033 

300.474 

300.656 -0. 1 46E-04 

15.000 

0.6933 

.4027 

.4022 

.4020 

300.714 

300.897 -0.21 9E-04 

20.000 

0.9314 

.4014 

.4009 

.4008 

300.956 

301.141 -0.293E-04 

25.000 

1.1718 

.4001 

.3996 

.3994 

301.201 

301.388 -0.366E-04 

30.000 

1.4146 

.3988 

.3983 

.3981 

301.448 

301.637 -0.440E-04 

35.000 

1.6598 

.3975 

.3970 

.3968 

301.697 

301.888 -0.51 3E-04 

40.000 

1.9075 

.3961 

.3957 

.3955 

301.949 

302.141 -0.587E-04 

45.000 

2.1577 

.3948 

.3943 

.3941 

302.204 

302.398 -0.661 E-04 

50.000 

2.4104 

.3934 

.3930 

.3928 

302.461 

302.656 -0.735E-04 

52.778 

2.5519 

.3927 

.3922 

.3920 

302.605 

302.801 -0.776E-04 


Final values 
Mole fractions: 

He « .2001 Ar = 0.0000 C02 = .3706 


CO= .0195 02 = .0098 N2 = .4001 
Gas Pressure (kPa): 101.325 

Reynolds number = 9.817 

Gas Heat Capacity (J/K-g) = 1.799 

Gas Density (g/cm A 3) = .0011720 

Gas Velocity (mm/s) = 39.378 

Effectivness factor = .99965 

Bulk Gas Diffusivity(cm A 2/s)= .22832 
Effective Diffusivity(cm A 2/s)= .053988 
Step size(mm) = .555556 


SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No.Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flowrate(liters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 
Computational loco parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): Length. Length(mm): 100.000 
Computation loop step size(mm): .5556 Display every 5.00 mm 


Initial values 

Reynolds number = 9.880 

Gas Heat Capacity (J/K-g) = 1.797 

Gas Density (g/cm A 3) = .0011817 

Gas Velocity (mm/s) = 39.063 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 

Effective Diffusivity(cm A 2/s)= .053407 

Step size(mm) = .555556 
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Distance 

%Conver 02gas 

02wall 02center 

Tgas 

Twall DPress 

(mm) 


<===(mMOL/L)= 

==> 

(Kelvin) 

(Kelvin) (kPa) 

0.000 

0.0000 

.4065 

.4060 

.4058 

300.000 

300.179 -0.000E+ 00 

5.000 

0.2239 

.4052 

.4048 

.4046 

300.236 

300.416 -0.730E-05 

10.000 

0.4574 

.4040 

.4035 

.4033 

300.474 

300.656 -0.1 46E-04 

15.000 

0.6933 

.4027 

.4022 

.4020 

300.714 

300.897 -0.21 9E-04 

20.000 

0.9314 

.4014 

.4009 

.4008 

300.956 

301.141 -0.293E-04 

25.000 

1.1718 

.4001 

.3996 

.3994 

301.201 

301.388 -0.366E-04 

30.000 

1.4146 

.3988 

.3983 

.3981 

301.448 

301.637 -0.440E-04 

35.000 

1 .6598 

.3975 

.3970 

.3968 

301.697 

301.888 -0.51 3E-04 

40.000 

1.9075 

.3961 

.3957 

.3955 

301.949 

302.141 -0.587E-04 

45.000 

2.1577 

.3948 

.3943 

.3941 

302.204 

302.398 -0.661 E-04 

50.000 

2.4104 

.3934 

.3930 

.3928 

302.461 

302.656 -0.735E-04 

55.000 

2.6657 

.3921 

.3916 

.3914 

302.720 

302.918 -0.809E-04 

60.000 

2.9236 

.3907 

.3902 

.3900 

302.983 

303.182 -0.883E-04 

65.000 

3.1842 

.3893 

.3888 

.3886 

303.247 

303.449 -0.957E-04 

70.000 

3.4475 

.3879 

.3874 

.3872 

303.515 

303.718 -0.1 03E-03 

75.000 

3.7136 

.3865 

.3860 

.3858 

303.785 

303.990 -0.11 IE-03 

80.000 

3.9825 

.3851 

.3846 

.3843 

304.058 

304.265 -0.1 1 8E-03 

85.000 

4.2543 

.3836 

.3831 

.3829 

304.334 

304.543 -0.1 26E-03 

90.000 

4.5290 

.3822 

.3817 

.3814 

304.613 

304.824 -0.1 33E-03 

95.000 

4.8066 

.3807 

.3802 

.3800 

304.895 

305.108 -0.141 E-03 

100.000 

5.0873 

.3792 

.3787 

.3785 

305.180 

305.395 -0. 1 48E-03 


Final values 
Mole fractions: 

He = .2001 Ar = 0.0000 C02 = .3712 


CO = .0190 02 = .0095 N2 = .4002 
Gas Pressure (kPa): 101.325 

Reynolds number = 9.753 

Gas Heat Capacity (J/K-g) = 1.800 

Gas Density (g/cm A 3) = .0011624 

Gas Velocity (mm/s) = 39.693 

Effectivness factor = .99961 

Bulk Gas Diffusivity(cm A 2/s)= .23125 

Effective Diffusivity(cm A 2/s)= .054566 

Step size(mm) = .555556 


SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No. Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas Composition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flowrate(liters/s): 2.500 Gas Temperature(K): 300.00 
Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 
Computational loop parameters: 

Output file (Full Profile/Summary): Full Profile 

Termination on (02 conversion/Length): Length. Length(mm): 100.000 
Computation loop step size(mm): 5.0000 Display every 5.00 mm 

Initial values 
Reynolds number = 

Gas Heat Capacity (J/K-g) = 

Gas Density (g/cm A 3) = 


98.802 

1.797 

.0011817 



Gas Velocity (mm/s) = 390.625 

Effectivness factor = .99969 

Bulk Gas Diffusivity(cm A 2/s)= .22540 
Effective Diffusivity(cm A 2/s)= .053407 

Step size(mm) = 5.000000 


Distance 

0) 

> 

c 

o 

O 

-0 

02gas 02wall 

02center 

Tgas 

Twall DPress 

(mm) 


<== 

=(mMOL/L)===> 

(Kelvin) 

(Kelvin) (kPa) 

0.000 

.0000 

.4065 

.4060 

.4058 

300.000 

300.179 -0.000E+00 

5.000 

.0156 

.4064 

.4059 

.4057 

300.024 

300.203 -0.730E-04 

10.000 

.0391 

.4062 

.4058 

.4056 

300.047 

300.227 -0. 1 46E-03 

15.000 

.0624 

.4061 

.4057 

.4055 

300.071 

300.250 -0.21 9E-03 

20.000 

.0855 

.4060 

.4055 

.4054 

300.095 

300.274 -0.292E-03 

25.000 

.1086 

.4059 

.4054 

.4052 

300.119 

300.298 -0.365E-03 

30.000 

.1317 

.4057 

.4053 

.4051 

300.142 

300.322 -0.438E-03 

35.000 

.1549 

.4056 

.4052 

.4050 

300.166 

300.346 -0.51 IE-03 

40.000 

.1781 

.4055 

.4050 

.4049 

300.190 

300.370 -0.584E-03 

45.000 

.2013 

.4054 

.4049 

.4047 

300.214 

300.394 -0.657E-03 

50.000 

.2245 

.4052 

.4048 

.4046 

300.238 

300.418 -0.730E-03 

55.000 

.2478 

.4051 

.4047 

.4045 

300.261 

300.442 -0.803E-03 

60.000 

.2711 

.4050 

.4045 

.4043 

300.285 

300.466 -0.877E-03 

65.000 

.2944 

.4049 

.4044 

.4042 

300.309 

300.490 -0.950E-03 

70.000 

.3177 

.4047 

.4043 

.4041 

300.333 

300.514 -0.1 02E-02 

75.000 

.3411 

.4046 

.4041 

.4040 

300.357 

300.538 -0.110E-02 

80.000 

.3644 

.4045 

.4040 

.4038 

300.381 

300.562 -0.117E-02 

85.000 

.3879 

.4043 

.4039 

.4037 

300.405 

300.587 -0.124E-02 

90.000 

.4113 

.4042 

.4038 

.4036 

300.429 

300.611 -0.1 32E-02 

95.000 

.4347 

.4041 

.4036 

.4034 

300.453 

300.635 -0.139E-02 

100.000 

.4582 

.4040 

.4035 

.4033 

300.477 

300.659 -0.1 46E-02 


Final values 
Mole fractions: 


He = .2000 Ar = 0.0000 C02 
CO = .0199 02 = .0100 N2 = 
Gas Pressure (kPa): 1 

Reynolds number = 

Gas Heat Capacity (J/K-g) = 
Gas Density (g/cm A 3) = 

Gas Velocity (mm/s) = 
Effectivness factor = 

Bulk Gas Diffusivity (cm A 2/s)= 
Effective Diffusivity (cm A 2/s)= 
Step size(mm) = 


*= .3701 
; .4000 
01 .324 
98.705 
1.798 
.0011800 
391.188 
.99969 
.22592 
.053513 
5.000000 


SELECTABLE PARAMETER SUMMARY 
Monolith Dimensions(mm): 

Support wall thickness: 1.00 Face dimension: 100.00 x 100.00 

Channel inner dimension: 4.00 No.Face channels: 20.00 x 20.00 
Active layer thickness: .25 % monol.volume open : 64.0 

Monolith inlet parameters: 

Gas ComDosition (mole fraction): C02: .3700 CO: .0200 02: .0100 

N2: .4000 He: .2000 Ar: 0.0000 
Gas Flowrate(liters/s): .250 Gas Temperature(K): 300.00 

Inlet Gas Pressure (kPa): 101.325 

Catalyst Properties: 

Catalyst Density (g/mm A 3): 0.500E-03 

Reaction rate constant at 298K(mm A 3 /gcat-s): 123.40 

Activation energy(J/mol): 39700.00 

Void-fraction as micropores: .24 Void-fraction as macropores: .48 
Avg. micropore radius (nm): .12E+02 Avg. macropore radius (nm): .50E+03 
Thermal Operation (adiabatic/isothermal): Adiabatic 
Computational loop parameters: 

Output file (Full Profile/Summary): Summary 

Termination on (02 conversion/Length): Length. Length(mm): 100.000 


< Q I 


Computation loop step size(mm): .5556 


Initial values 
Reynolds number = 

Gas Heat Capacity (J/K-g) = 

Gas Density (g/cm A 3) = 

Gas Velocity (mm/s) = 

Effectivness factor = 

Bulk Gas Diffusivity(cm A 2/s)= 

Effective Diffusivity(cm A 2/s)= 

Step size(mm) = 

Distance %Conver 02gas 02wall 02center Tgas Twall DPress 
(mm) <===(mMOL/L)===> (Kelvin) (Kelvin) (kPa) 

0.000 0.0000 .4065 .4060 .4058 300.000 300.179 -0.000E+00 

100.008 5.0878 .3792 .3787 .3785 305.180 305.395 -0.148E-03 

Final values 
Mole fractions: 

He = .2001 Ar = 0.0000 C02 = .3712 
CO = .0190 02 = .0095 N2 = .4002 
Gas Pressure (kPa): 101.325 

Reynolds number = 9.753 

Gas Heat Capacity (J/K-g) = 1.800 

Gas Density (g/cm A 3) = .001 1 624 

Gas Velocity (mm/s) = 39.693 

Effectivness factor = .99961 

Bulk Gas Diffusivity(cm A 2/s)= .23125 

Effective Diffusivity(cm A 2/s)= .054566 

Step size(mm) = .555600 


File CATONE. See if you can pick out some of the parameters listed below. The program 
listing will show the structure of the parameter file. 

100.0 

100.0 

.5000 

4.000 

4.000 
.2500 
.2000 

1.0000E-07 

.3700 

2.0000E-02 

1.0000E-02 

.4000 

1.000 


2.500 

100.0 

3.9700E+04 

123.4 

300.0 
.2500 
9 

.5556 

5.0000E-04 

.2400 

.4800 

12.00 

500.0 


9.880 

1.797 

.0011817 

39.063 

.99969 

.22540 

.053407 

.555600 
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Appendix 3; Menu Description 


LASCAT Main Menu 

The LASCAT Main Menu section provides detailed instructions on the use of the main 
menu. It is assumed that the LASCAT Tutorial has been reviewed. The LASCAT Main Menu is 
first encountered upon program execution and has five options available. Each option is detailed 
below. To select an option type the desired choice (1- 5) and then RETURN (key). A number 
other than 1-5 will cause an input requirement message to appear and another chance to enter 1-5 is 
made available. 

Option 1 - Read in new operating parameters: 

If you have a file containing a selectable parameter set that you want to use, choose this 
option. This file would have been previously generated by saving the parameter set using 
option 15 of the parameter modification menu. You will first be queried whether the file is 
New (N) or Existing (E). the file is Existing (E). If you indicate that the file is existing, 
you will be queried for the file name. Some systems are sensitive to upper and lowercase 
differences in file names. After specifying the filename, you are returned to the main 
menu. If you indicate that the file is new, default file parameters are assigned. 

Option 2 - Show current operating parameters: 

This option lists the current selectable parameter set. Use this option to verify the proper 
selectable parameter set prior to execution of the program. After parameter set review, you 
arc returned to the main menu. 

Option 3 - Change operating parameters: 

This option leads to the parameter modification menu described below. The options similar 
to main menu options 1-5 are available in the parameter modification menu. You can 
return to the main menu from the parameter modification menu. 

Option 4 - Run program: 

The Run program option executes the computational portion of the program, using the 
current selectable parameter set. Results obtained are output to the screen for review and to 
the file LDATA for further processing or review. Any additional "runs" in the session will 
be tacked on to the end of file LDATA. (Macintosh only: If the program is terminated and 
rerun the old LDATA file will be erased and a new LDATA file started . If you desire to 
keep the information in contained in LDATA, rename LDATA prior to running LASCAT 
again.) After computations are completed or a termination condition is met, the program 
returns to the main menu. 

Option 5 - Exit proeram: 

The program is exited. If using the Macintosh remember the caution concerning LDATA, 
as noted above in the Option 4 section. 


LASCAT Parameter Modification Menu 

The LASCAT parameter modification menu section provides detailed instructions on the 
use of the parameter modification menu. It is assumed that the LASCAT Tutorial has been 
reviewed. The LASCAT parameter modification menu is encountered when option 3 of the main 
menu is chosen and has eighteen options available. Each option is detailed below. To select an 
option type the desired choice (1-18) and then RETURN. A number other than 1-18 will cause an 
input requirement message to appear and another chance to enter 1-18 is made available. The 
default file parameter set is the parameter set is instituted at the start of the program. 
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Option 1 - Change monolith physical dimensions: 

A number of physical dimensions (see diagram below) are queried for in this section. A 
description of each dimension and the associated restrictions are listed. Values that 
represent a valid physical situation and meet the listed restrictions should not cause any 
problems. A number of constraint for given below. The program assumes monolith 
channels are square. Remember that the system under study can be scaled up or down. 
For instance, a monolith with an inlet gas flow rate of 1 liters/sec and a face (H x V) of 
100 x 100 mm can be scaled up to an inlet gas flow rate of 4 liters/sec and a face of 200 x 
200mm. The other dimensions (H,W, and A) are held the same. 

Monolith horizontal face dimension (mm), H, 10 < H < 1000 
Monolith vertical face dimension (mm), V, 10 < V £ 1000 
Monolith channel wall thickness (mm), T, T <H, T <V 
Monolith channel opening dimension (mm) , W, W < H, W < V 
Monolith active layer thickness (mm) , A , A < W/2 

Default file dimension values(mm) are H=V= 100.0, T=0.50, W=4.00, and A =0.25 
Monoliths with square channels but with circular, oval, etc., shaped faeces can be 
modelled. For these nonstandard shapes, determine the dimensions ( H x V ) of a square 
faced monolith that has an equivalent facial area and use the calculated H and V in the 
program. 

Option 2 - Change inlet gas composition: 

The program allows for the laser operating gas to contain He, Ar, CO 2 , O 2 , CO, and N 2 
gases. The species mole fraction is defined as the mole percent divided by 100. For 
example, 10 mole% N 2 corresponds to a mole fraction of 0.10. Each mole fraction must 
be between 0 and 1 . Mole fractions for He, Ar, CO 2 , O 2 are quieried for when this option 
is chosen. Stoichiometric ratios of C> 2 :CO (1:2) are assumed; the input mole fraction of O 2 
is used to calculate the CO mole fraction. The mole fraction of O 2 is limited to 0.03 
(3%).The sum of the mole fractions must add up to 1; N 2 mole fraction is calculated using 
the mole fractions of the other species. The default file values for mole fractions are 
He=0.20, Ar=0.0, CO 2 =0.37, O 2 =0.01, CO=0.02, and N 2 =0.40. 

Option 3 - Change inlet gas volumetric flow rate: 

Inlet gas volumetric flow rate (liters/sec) is queried for in this option. The inlet flow rate is 
assumed to be at the temperature specified in option 4. Although there is no upper limit on 
the inlet gas flow rate, if the flow rate results in a Reynolds (Re) number > 2300, the 
computational portion of the program will be terminated. The default file value is 0.25 
liters/sec. 

Option 4 - Change inlet gas temperature: 

Inlet gas temperature (K) is quieried for in this option. The acceptable range of temperature 
is 200 to 700 K. Gas temperature is updated along the length of the monolith. If the gas 
temperature reaches 700 K along the length of the monolith the computational portion of the 
program will be terminated. The default file inlet gas temperature is 300 K. 

Option 5 - Change inlet gas pressure: 

Inlet gas pressure (kPa) is quieried for in this option. The acceptable range for inlet gas 
pressure is 0.1 to 350 kPa. The pressure is updated along the length of the monolith. If 
the gas pressure exceeds 350 kPa or drops below O.lkPa along the length of the monolith 
the computational portion of the program will be terminated. The default file value is 
101.325 kPa ( latm ). 

Option 6 - Change catalyst activation energy: 

Catalyst activation energy (J/mol) is quieried for in this option. The acceptable range for 
catalyst activation energy is between 6 and 300,000 J/mol. The default file value is 39700 
J/mol . 

Option 7 - Change catalyst reaction rate constant: 
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Catalyst reaction rate constant (mm 3 /gcat-sec) [referenced to 298 K] is quieried for. The 
catalyst activation energy must be greater than zero. The default file value is 123.4 
mm 3 /gcat-sec. 

Option 8 - Change catalyst active layer density: 

Catalyst active layer density (g/mm 3 ) is quieried for in this option. The catalyst active layer 
density must be greater than zero. The default file value is 5e-04 g/mm 3 . 

Option 9 - Change macro/micro pore radius & void-fraction: 

Four quantities, macropore and micropore radius and macropore and micropore void- 
fractions, are queried for in this option. 

Define macropore radius (nm) = RMACRO Default value= 500 nm 
Define micropore radius (nm) = RMICRO Default value= 12 nm 
Define macropore void-fraction = EMACRO Default value= 0.48 
Define micropore void-fraction = EMICRO Default value= 0.24 
Acceptable ranges: RMICRO < RMACRO, EMACRO + EMICRO < 1.0 

Option 10 - Change thermal operation (Adiabatic/Isothermal 1: 

Type of monolith thermal operation (adiabatic/isothermal ) is quieried for in this option. For 
isothermal operation the gas and the monolith section are maintained at (or very close to)the 
gas inlet temperature specified in option 4. Default file operation mode is adiabatic. 

Option 1 1 - Change outp ut profile (Full/ Summary ): 

Type of output profile is quieried for in this option. Only the initial and final values 
(temperature, pressure, and O2 concentration) are output in the summary option. In the full 
option, the values are specified for each display interval as specified in option 13. Default 
file value is Full. See LASCAT Tutorial for an example. 

Option 12 - Change termination variable (Con version/Length ): 

The termination condition (O2 conversion or monolith length) is queried for in this option. 
If O2 conversion is chosen, the termination percent conversion ( 0-100 %) is required to be 
specified The default O2 conversion is 2.5 % If the monolith length option is chosen, the 
monolith length (5-4000mm) is required to be specified. The default monolith length is 
5mm. See LASCAT Tutorial for an example. 

Option 13 - Change computation display interval: 

Computaion display interval (CDI) is quieried for in this option. The display interval 
specifies the reporting interval for temperature, pressure, and O2 concentration. The 
specified display interval may be over overriden if too many computational loop iterations 
(>250) are required to meet the specified display interval. The default file value is 5mm. 

Option 14 - Show current operating parameters: 

This option lists the current selectable parameter set. Use this option to verify the proper 
selectable parameter set prior to execution of the program. After parameter set review, you 
are returned to the parameter modification menu. This option is similar to option 2 of the 
main menu. 

Option 1 5 - Save current operating pararaeterc; 

This option allows the user to save the current operating parameter set to file. The file 
name and the file status (new or existing) are required to be specified. The main menu's 
option 1 allows parameter set files to be installed as the current selectable parameter set. 

Option 1 6 - Run program: 

The "Run program" option executes the computational portion of the program, using the 
current selectable parameter set. Results obtained are output to the screen for review and to 
the file LDATA for further processing or review. Any additional "runs" in this session will 
be tacked on to the end of file LDATA. (Macintosh only: If the program is terminated and 
rerun the old LDATA file will be erased and a new LDATA file started . If you desire to 
keep the information in contained in LDATA, rename LDATA prior to running LASCAT 
again.) After computations are completed or a termination condition is met, the program 
returns to the main menu. This option is identical to option 4 of the main menu. 

Option 17 - Return to main menu: 
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Returns to the main menu. See main menu seciton above. 

Option 1 8 - Exit program; 

The program is exited. If using the Macintosh remember the caution concerning LDATA, 
as noted above. 
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1 £************************ program lascat************************ 

2 C Written by 

3 C Keith Guinn, Seth Goldblum Dr. Richard Herz,and Ed Noskowski 

4 C Department of AMES / Chemical Engineering B-010 

5 C University of California, San Diego 

6 C La Jolla, CA 92093 

7 C 

8 C Date revised-February 1989 

9 C 

10 C Funded by NASA's Langley Research Center GRANT No. NAG-1-823 

11 C 

12 C Fortran Software Compatibility 

13 C MICROSOFT FORTRAN Version 2.2 (MACINTOSH PLUS,SE, & II) 

14 C ABSOFT MacFORTRAN / 0 2 0 Version 2.3 (MACINTOSH II) 

15 C This program is FORTRAN 77 compatible. 

16 C 

17 C Editing program used: QUED Version 1.53 (MACINTOSH) 


18 

c 




19 

c 

Program listing outline 


line numbers 

20 

c 

1 . Program information and table of 

contents 

2 -35 

21 

c 

1. Indentifiers 


36 -153 

22 

c 

2 . Program Algorithm/Restrictions 


154 -297 

23 

c 

3 . Program LASCAT 


298 -1450 

24 

c 

A. Variable declaration 


299 -322 

25 

c 

B . Introduction/purpose 


323 -357 

26 

c 

C. Default/read/write parameters 


358 -457 

27 

c 

D . Main-menu, Sub-menu 


457 -551 

28 

c 

E. Parameter input section, write 

parameters 

552 -981 

29 

c 

F. Initialize physical properties 


982 -1038 

30 

c 

G. Step size check loop 


1039 -1116 

31 

c 

H. Main computational loop 


1162 -1409 

32 

c 

4. Subroutines: RK4,FUNC, LINES, HEADER 

1451 -1906 

33 

c 

TCON, PROP, VIS, FXN1 



34 

c 

FXN2 , FXN3 , FXN4 



35 

c 





36 c******************** IDENTIFIERS ★******************************* 

37 C A-REACTION RATE CONSTANT OF ACTIVE CATALYST MATERIAL (mm~3/c-cat s) 

38 C AMW-AVERAGE MOLECULAR WEIGHT OF GAS MIXTURE (g/mole) 

39 C AREF -VOLUMETRIC PUMPING RATE OF ACTIVE CATALYST AT 298K(mm A 3/g-cat S) 

40 C ASK-QUERY DIRECTION NUMBER 

41 C AX (i, j) -i, j=l->6** CROSS TERM MATRIX FOR THERMAL CONDUCTIVITY 

42 C BEEP-BEE? CHARACTER 

43 C CDICOMPUTATION DISPLAY INTERVAL (mm) 

44 C CGX-MONOLITH CHANNEL 02 GAS CONC. AT POSITION X (mol/mm~3) 

45 C CGXTMP-MONOLITH CHAN. 02 GAS CONC. AT POSITION X-DX (mol /mm"' 3) 

46 C THIS IS A HOLDING CONSTANT TO DETERMINE DCGX 

47 C CGXZ -MONOLITH CHAN. 02 GAS CONC. AT POSITION X=0 (mol /mm'' 3) 

48 C CH, CHI, CH2 - CHARACTER* 1 Dummy Variable 

49 C CHFLO-MONOLITH CHANNEL VOLUMETRIC FLOWRATE (mm~3/s) 

50 C CHGVEL-MONOLITH CHANNEL GAS VELOCITY (mm/s) 

51 C CKHD-MONOLITH CHANNEL HYDRAULIC DIAMETER (mm) 

52 C CHSTV -MONOLITH CHANNEL SURFACE AREA TO VOLUME (1/mm) 

53 C CHWP-MONOLITH CHANNEL WETTED PERIMETER (mm) 

54 C CMCWT-CHARACTERISTIC MCWT (mm) 

55 C CONST1 , CONST2 , CONST3 , CONST4 , CONST5-RESULT HOLDING CONSTANTS 

56 C CONST5, CONST6, CONST7, CONST8, CONST9-RESULT HOLDING CONSTANTS 

57 C CONSTA, CONSTB, CONSTD, CONSTE, CONSTF, CONSTC-RESULT HOLDING CONSTS 


Appendix 4 - program listing 
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58 C CONV-CONVERSION OF OXYGEN IN MONOLITH CHANNEL (%) 

59 C CONVX-CONVERSION OF OXYGEN AS FUNCTION OF X (%) 

60 C CONVXZ-CONVERSION OF OXYGEN AT X=0 (%) 

61 C CP (i),i=l, 6-HEAT CAPACITY OF SPECIES i AT TEMP. TGX (J/K-g) 

62 C CPGAS-HEAT CAPACITY OF BULK GAS MIXTURE (J/K-g) 

63 C CPCAT-HEAT CAPACITY OF CATALYST WASHCOAT (J/K-g) 

64 C CPSUP-HEAT CAPACITY OF MONOLITH SUPPORT (J/K-g) 

65 C CSX-MONOLITH CHAN. WALL 02 GAS CONC AT POSIT. X (mol/mm A 3) 

66 C CSXZ-MONOLITH CHAN. WALL 02 GAS CONC AT POSIT X=0 (mol/mm A 3) 

67 C D(i) , i=l- 6-BINARY DIFF COEFF FOR 02 in i 298K,latm (mm A 2/s) C 

68 C DAB-DIFFUSIVITY OF 02 IN BULK GAS (mm A 2/s) 

69 C DABEFF-EFFECTIVE DIFFUSIVITY OF 02 IN WASHCOAT (mm A 2/s) 

70 C DCGX-CHANGE IN CGX FOR A CHANGE DX(mol/mm A 3) 

71 C DENGAS -DENSITY OF BULK GAS MIXTURE (g/mm A 3) 

72 C DENGSR-REFERENCE DENSITY OF BULK GAS MIXTURE (g/mm A 3) 

73 C DENCAT-DENSITY OF CATALYST WASHCOAT (g/mm A 3) 

74 C DENSUP-DENSITY OF MONOLITH SUPPORT (g/mm A 3) 

75 C DP-PRESS CHANGE ASSOC. WITH A CHANGE IN CHANNEL POSITION DX(atm) 

76 C DFGX- CHANGE IN FGX FOR A CHANGE DX (MOL/s) 

77 C DKMI CR-KNUDSEN DIFFUSION COEFFICIENT IN MICROPORES (mm A 2/s) 

78 C DKMACR-KNUDSEN DIFFUSION COEFFICIENT IN MACROPORES (mm A 2/s) 

79 C DMICRO-EFFECTIVE DIFFUSION COEFFICIENT IN MICROPORES (mm A 2/s) 

80 C DMACRO-EFFECTIVE DIFFUSION COEFFICIENT IN MACROPORES (mm A 2/s) 

81 C DTGX-CHANGE IN TGX FOR A CHANGE DX(mol/mm A 3) 

82 C DX-INCREMENTAL CHANGE IN MONOLITH CHANNEL POSITION X (mm) 

83 C DXC HK- COUNTS # OF TIMES THROUGH THE LOOP DURING STEP SIZE CHECK 

84 C EMICRO-FRACTION OF TOTAL VOIDS IN PELLETS AS MICROPORES 

85 C EMACRO-FRACTION OF TOTAL VOIDS IN PELLETS AS MACROPORES 

86 C ENGACT-CATALYST REACTION ACTIVATION ENERGY ( J/mol ) 

87 C FZERO - INLET 02 MOLAR FLOWRATE (MOL/s) 

88 C FACTOR-PROFILE STEP SIZE FACTOR (integer) 

89 C FGX - MOLAR FLOWRATE OF 02 AT POSITION X (MOL/s) 

90 C FGXEND - OUTLET 02 MOLAR FLOWRATE (MOL/s) 

91 C FILNAM - NAME OF DATAFILE CONTAINING DEFAULT PARAMETERS 

92 C FIRST - LOGICAL* 1 IS TRUE IF FIRST TIME THROUGH LOOP, AND 
S3 C FALSE IF ANYTHING ELSE 

94 C H-HEAT TRANSFER COEFFICIENT FROM GAS TO WALL ( J/mm''2-S-K) 

95 C HEATRX- NEGATIVE CHANGE IN ENTHALPY FOR 02+2CO->CO2 REACTION (J/mol) 

96 C I, J-INTEGER COUNTER 

97 C JUMP-INTEGER- VALUE DETERMINES IF DX CAN BE FURTHER REDUCED 

98 C K-TKSRMAL GAS CONDUCTIVITY (J/mm-S-K) 

99 C KC-MASS TRANSFER COEFFICIENT (mm/ s) 

ICC C KREF -REFERENCE THERMAL GAS CONDUCTIVITY ( J/mm A 2-s-K) 

101 C LAMBMX - THERMAL CONDUCTICITY OF THE BULK GAS MIXTURE (mW/cm K) 

1C2 C MAIN - 300LEAN EXPRESSION. TRUE WHEN USER IS IN MAIN MENU 

103 C MCH-MONOLITH CHANNEL HEIGHT (mm) 

104 C MCST-MONOLITH CHANNEL SUPPORT WALL THICKNESS (mm) 

105 C MCW-MONOLITH CHANNEL WIDTH (mm) 

106 C MCWT -MONOLITH CHANNEL WASHCOAT THICKNESS (mm) 

107 C MFDH-MONOLITH HORIZONTAL FACE DIMENSION (mm) 

108 C MFDV-MONOLITH VERTICAL FACE DIMENSION (mm) 

109 C MOLFLO - MOLAR FLOWRATE OF BULK GAS THROUGH MONOLITH (M/s) 

110 C MONOL-COMPUTATION LOOP COMPLETION MONOLITH LENGTH (mm) 

111 C MVFLO-MONOLITH VOLUMETRIC GAS FLOWRATE (mm A 3/ s) 

112 C MVFLZ -INLET MONOLITH VOLUMETRIC GAS FLOWRATE (mm A 3/s) 

113 C MW (i) -i->l, 6** MOLECULAR WEIGHT OF COMPONENT GASES (gms /mole) 

114 C N-EFFECTIVESS FACTOR FOR CATALYST WASHCOAL (NONDIM) 



Appendix 4 - program listing 


ORIGINAL PAuE I** 

OF POOR QUALITY 


115 C NOCHM-NUMBER OF CHANNELS IN MONOLITH 

116 C NOCHV-NUMBER OF CHANNELS VERTICALLY IN MONOLITH 

117 C NOCHH-NUMBER OF CHANNELS HORIZONTALLY IN MONOLITH 

118 C 02ZER0-M0N0LITH INLET OXYGEN CONCENTRATION (%02) 

119 C ONCE-INTEGER-DETERMINES IS THIELE STATEMENT HAS BEEN WRITTEN 

120 C P-DO LOOP VARIABLE 

121 C PZERO-INLET PRESSURE OF GAS MIXTURE (atm) 

122 C PRESS-PRESSURE OF GAS MIXTURE IN THE CHANNEL (atm) 

123 C Q-HOLDING CONSTANT 

12- 9 C R-GAS CONSTANT (J/mol-K) 

125 C RE - Reynold’s number (dimensionless) 

126 C RMICRO- AVERAGE MICROPORE RADIUS (cm) 

127 C RMACRO- AVERAGE MACROPORE RADIUS (cm) 

128 C S-TAB CHARACTER 

129 C SS-CARRIAGE RETURM CHARACTER 

130 C SSS-A OR I FOR ADIABATIC OR ISOTHERMAL-CHARACTER 

131 C SSS1-P OR S FOR PROFILE OR SUMMARY-CHARACTER 

132 C SSS2-0 OR L FOR 02 CONVERSION OR MONOLITH LENGTH- CHARACTER 

133 C T-GAS TEMPERATURE FOR VISCOSITY AND THERMAL CONDUCTIVITY 

13- 3 C TCON(i, j)-i->l,6**, j->2, 7***, THERMAL COND. OF COMPONENT GASES 

135 C AT VARIOUS TEMPERATURES (mW/cm-K) 

136 C TCONT (i ) -i->l , 6**THERM COND. OF COMPONENT GASES, @TEMP T (mW/cm*K) 

137 C TGX -MONOLITH CHANNEL BULK GAS TEMPERATURE AT POSIT. X (K) 

138 C THIELE-SLAB THIELE MODULUS FOR CATALYST WASHCOAT (NONDIM) 

139 C TREF-PHYSICAL PROPERTIES REF TEMP=298K (K) 

1-90 C TSX -MONOLITH CHANNEL WALL TEMPERATURE AT POSIT. X (K) 

141 C TSXZ -MONOLITH CHANNEL WALL TEMPERATURE AT POSIT. X=0 (K) 

142 C TZERO-MONOLITH INLET GAS TEMPERATURE (K) 

143 C VISC (i, j)-i->l, 6**, j->2, 7***, VISCOSITY OF COMPONENT GASES 

144 C AT VARIOUS TEMPERATURES (centipoises) 

145 C VISCT (i) -i->l, 6** VISCOS OF COMP GASES @ TEMP T (centipoises) 

146 C VISCTC-COMPOSITE GAS VISCOSITY (centipoises) 

147 C X -MONOLITH CHANNEL POSITION (mm) 

148 C XZ-MONOLITH CHANNEL POSITION AT X=0 (mm) 

149 C Y (i) -i->l, 6** MOLE FRACTION OF COMPONENT GASES (gms /mole) 

150 C Y10, Y20, Y30, Y40, Y 50, Y 60-HOLDING CONSTANTS FOR Y(i)'s 

151 C **l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2 

152 C ***2=200, 3=300, 4=400, 5=500, 6=600, 7=700K 

153 C 

154 c*********** PROGRAM ALGORITHM/RESTRICTIONS ******************** 

155 C -Declare variables 

156 C -Open output data file "LDATA" using unit 25 

157 C -Display introduction on screen 

158 C -Read in values of default operating parameters 

159 C -Compute some geometric and record starting values using default 

160 C parameters 

161 C -Display main menu, query for desired option (1-5) 

162 C 1-Read in new operating parameters 

163 C -Query for file from which to read from 

164 C -Read in values from file 

165 C -Return to main menu 

166 C 2 -Show current operating parameters 

167 C -Display operating parameter summary table, this includes some 

168 C calculated values 

169 C -Return to main menu 

170 C 3-Change operating parameters 

171 C -Drop to* submenu, query for desired option (1-18), if not 1-18 
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172 C 

173 C 

174 C 

175 C 

176 C 

177 C 

178 C 

179 C 

180 C 

181 C 

182 C 

183 C 

184 C 

185 C 

186 C 

187 C 

188 C 

189 C 

190 C 

191 C 

192 C 

193 C 

194 C 

195 C 

196 C 

197 C 

198 C 

199 C 

200 C 

201 C 

202 C 

203 C 

204 C 

205 C 

206 C 

207 C 

208 C 

209 C 

210 C 

211 C 

212 C 

213 C 
•214 C 

215 C 

216 C 

217 C 

218 C 

219 C 

220 C 

221 C 

222 C 

223 C 

224 C 

225 C 

226 C 

227 C 

228 C 


query again. After each section is completed return to 
the submenu. 

-1 Change monolith physical dimensions (mm) 

-Input monolith face dimensions (vertical & horizontal) 
range: >=10mm, <=10 00mm 

-Input monolith channel wall thickness, inner channel 
dimension, active layer thickness 
range: support wallthickness 

< either vert, or horiz face 
range: inner channel dimension 

< either vert or horiz face 

range: active layer thickness < 0.5 * ( support wall) 
-Recompute some geometric parameters using new dimensions 
-2 Change Inlet Gas Composition (mole fractions) 

-Input mole fractions of He, Ar, C02, and 02 
-Mole fraction of CO is set to twice that of 02 
-Mole fraction of N2 is computed to make mole fraction 
sum equal to one 

range: each mole fraction >=0 and <=1, 
sum mole fractions =1 

-3 Change inlet gas volumetric flowrate (1/s) 

-Input inlet volumetric flowrate 
range: volumetric flowrate >0 
-4 Change inlet gas temperature (Kelvin) 

-Input inlet gas temperature 

range: temperature >=200K and <=700K 
-5 Change inlet gas pressure (kPa) 

-Input inlet gas pressure (kPa) 
range: pressure 0.1-350 kPa 
-6 Change catalyst activation energy (J/mol) 

-Input caltalyst activation energy (J/mol) 
range: activation energy >0 and <300000 
-7 Change catalyst reaction rate constant (mm A 3/gcat-s) 

-Input reaction rate constant 
range: reaction rate >0 

-8 Change catalyst active layer density (g/mm A 3) 

-Input density 

range : density >0 

-9 Change macro/micro pore radius & void-fraction 
-Input macro/micro pore radius 

range: radius >0 , micro<macro 
-Input macro/micro void-fractions 

range: sum <=1 , void-fractions>=0 
-10 Change thermal operation (Adiabatic/Isothermal) 

-Input option A/I or a/i 
-11 Change output profile (Full/ Summary) 

-Input option P/S or p/s 

-12 Change termination variable (Conversion/Length) 

-Input termination option 02 conversion (0) or monolith 
monolith lenght (L) . If (0) is chosen, specific 
termination conversion (%02) . If (L) is chosen specify 
monolith length (mm). 02% conversion range >0,<100%. 

Monolith length range >5, <4000mm. 

-13 Change computation display interval (mm) 

-Input computation display interval (mm) 

range: computation display interval <=100mm 
-14 Show current operating parameters 
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229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 
269 

i I U 

L. ! 1 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 


C -Displays current operating parameters 

C -15 Save current operating parameters 

C -Query for name of file which current operatings parameters 

C will be saved. Unit 24 is used. 

C -16 Run program 

C -Executes computational loop until desired termination 

C condition is reached. See below for details. 

C -17 Return to main menu 

C -See above for options in main menu 

C -18 Exit program 

C -Terminates program 

C 4 -Run program 

C -Executes computational loop until desired termination condition 

C is reached. See below for details. 

C 5-Exit program 
C 

C -Run program details: 

C -Computation of average molecular weight and gas density 
C -Computation of parameters needed for thermal conductivity and 
C viscosity calculations 

C -computation of channel volumetric flowrate and channel gas velocity 
C -initialize some physical properties parameters 
C -step size computational loop ( 2 loop passes) 

C -compute gas density, channel flow and velocity 
C -compute bulk diffusivity andthermal conductivity 

C -compute viscosity of gas and heat capacity of gas 

C -compute transport parameters and rxn rate 
C -compute thiele modulus and effectiveness factor 
C -compute Reynolds number and conversion 

C -perform Runge-Kutta integration 

C -update variables 

C -return to beginning of step size computational loop 
C -compute step size to produce approx .04% change per loop iteration 
C -initialize some physical properties parameters 

C -main computational loop (continue until termination criteria met) 

C -check if termination criteria are met, if so exit loop 
C -compute gas density, channel flow and velocity 
C -compute bulk diffusivity andthermal conductivity 

C -compute viscosity of gas and heat capacity of gas 
C -compute transport parameters and rxn rate 

C -compute thiele modulus and effectiveness factor 

C -compute Reynolds number and conversion 

C -print (to screen and file) header if first time through 

C -print (to screen and file) parameters at desired interval 

C -perform Runge-Kutta integration 

C -update variables 

C -step size check section, if conversion per step is < 0.08%, then 
C reduce step size (/5) and increase display factor (*5) 

C -parameter value boundary check, if out of bounds, exit loop 
C -return to beginning of main computational loop 

C -final computed values output section 
C -reset some important variable values 
C -return to main menu 
C -end of program 
C -subroutines 

C -RK4, fourth order Runge-Kutta 

C -FUNC, support for RK4 


\ 
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286 C -LINES, writes a specified no. of lines 

287 C -HEADER, writes a portion of the header to screen and file 

288 C -TCON, computes thermal cond. of each species and a specifed temp 

289 C -PROP, computed viscosity and thermal conductivity of gas 

290 C -VIS, computes the viscosity of each species at a specified temp 

291 C -OUTPUT, write a portion of the header to screen and file 

292 C -FXN1 , computes effective diffusivity of gas mixture 

293 C -FXN2, computes heat capacity of gas mixture 

294 C -FXN3, computes diffusivity of 02 in mixture 

295 C -FXN4, computes average molecular weight 

296 C 

297 

298 C PROGRAM LASCAT 

300 COMMON /BLK1/ KC, CHSTV,MCH,MCW, CSX, CGX, CHFLO, CPGAS 

301 COMMON /BLK2/ VISCTC,CHGVEL,CHHD,H,TSX,DENGAS 

302 INTEGER FACTOR, ONCE 

303 REAL* 8 MFDH,MCW,MCH,MCWT, CMCWT,MCST, CHWP, CHSTV, CHHD 

304 REAL* 8 TZERO, 02ZER0, MVFLO, CONV, CHFLO, CHGVEL, MONOL, SUM 

305 REAL* 8 CGX,CSX,TGX,TSX,TREF,X,DX,DCGX,DTGX,AREF,A,CDI 

306 REAL* 8 KC, DENGAS, DENGSR, DENCAT, DENSUP, CPGAS, CPCAT, CPSUP 

307 REAL* 8 N, HEATRX, ENGACT, R, THIELE, DAB, DABR, DABEFF, LAMBMX 

308 REAL* 8 H,K,KREF,CONSTl,CONST2,CONST3,P,Q,CONVX,ASK 

309 REAL*8 CONSTD,CONSTE,CONSTF,XO (3) ,XNEW(3) 

310 REAL*8 EMICRO, EMACRO, RMICRO, RMACRO, MFDV 

311 REAL* 8 DKMICR, DKMACR, DMICRO, DMACRO 

312 CHARACTER* 1 SS, SSS, SSS1, SSS2, BEEP, CH, CHI, CH2 

313 CHARACTER* 2 0 FILNAM, S 

314 REAL* 8 T, Y (6) , TCONDC, VISCTC, Y10, Y20, Y30, Y40, Y50, Y60 

315 REAL* 8 MW ( 6 ) , MOLFLO, MVFLZ , PZERO 

316 REAL*8 VISC (6, 2 : 7) , AX (6, 6) 

317 REAL*8 VISCT (6) , TCONT (6) , AMW, PRESS, D (6) , DP 

318 REAL* 8 CONST5, CONST4, CONST6, CONST7, CONST8, CP (6) 

319 REAL* 8 FGX, DFGX, FZERO, FGXEND, NOCHV, NOCHH, NOCHM 

320 REAL* 8 CONVXZ,XZ,CGXZ,CSXZ,TSXZ,RE 

321 INTEGER I, J, CONSTA, CONSTB, CONSTC, DXCHK,JUMP 

322 LOGICAL*! FIRST, MAIN 

323 c****** Introduction******** **********★★*****★■*■**★★■*•■* **★★**★★*★★** *** 

324 C Open Data file which will be used for output. 

325 OPEN ( UNIT=25, FILE='LDATA' , STATUS='NEW' ) 

326 WRITE (*, *) ' ****** PROGRAM LASCAT *★★★★★★» 

327 WRITE (*,*) 

328 WRITE (*, * ) 1 The purpose of this program is to calculate 1 , 

329 2 1 the gas concentration and 1 

330 WRITE (*, *) * temperature profiles of a monolith catalyst section 1 , 

331 2 1 of a C02 laser. The 1 

532 WRITE (* , * ) 1 C02 decomposes when the laser is pulsed. The CO ', 

333 2 1 and 02 produced as a result 1 

334 WRITE (*,*) ' of pulsing are detrimental to the efficient opera', 

335 2'tion of the laser. The' 

336 WRITE (*, *) 'recombination reaction is CO +1/2 02 -> C02. ', 

337 2 'This program provides the' 

338 WRITE (*,*) 'means to model the performance of a monolith 1 , 

339 2 'catalyst section under various' 

340 WRITE (*,*) 'gas compositions, temperatures, catalyst ', 

341 2 'activities gas flowrates, oxygen ' 

342 WRITE (*,*)' conversion, monolith face and length dimensions. ', 
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343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 


386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 


2' Results can indicate if ' 

WRITE (*,*)' constraints such as conversion, maximum gas temper', 
2'ature, monolith weight ' 

WRITE (*,*) 'are satisfied and how the system parameters may be', 
2' altered to meet these' 

WRITE (*,*)' constraints. Parameters and options may be altered', 
2' to tailor the monolith' 

WRITE (*,*) 'design. Default values can also be used as a ' , 

2 'starting point for the design' 

WRITE (*,*) 'process. A review of the parameters and options', 

2' chosen may be made prior' 

WRITE (*,*) 'to execution of the computational portion of the', 

2 ' program. ' 

CALL LINES (5) 

WRITE (*,*) ' (HIT RETURN TO CONTINUE) ' 

READ (*,621) S 

621 FORMAT (A) 

622 FORMAT (A1 ) 

BEEP= CHAR (7) 

C FIRST TIME THROUGH, SKIP OVER QUERY FOR READING PARAMETER FILE 

GO TO 2119 

C Default values for parameters 

C First see if operating parameters are already available in a data 

C file. 

9000 WRITE ( * , * ) 

WRITE (*,*) 'Read in data from existing parameter file ? (Y/N) ' 

READ (*,621) CH 

IF ( (CH.EQ. 'Y' ) .OR. (CH.EQ. 'y' ) ) THEN 

WRITE (* , * ) 

WRITE (*, *) 

WRITE (*,*) 'What is the name of the existing parameter file?' 

WRITE (*,*) 'to be opened?' 

READ (*, 621) FILNAM 

OPEN (UNIT=23, FILE=FILNAM, STATUS= ' OLD ' ) 

C If existing data file exists, read in data. 

READ (23, *)MFDH 
READ (23, *)MFDV 
READ (23,*) MOST 
READ (23, *)MCW 
READ (23, *)MCH 
READ (23, * ) MCWT 
READ (23, *) Y (1) 

READ (23, *)Y (2) 

READ (23, *) Y (3) 

READ (23, *) Y (4) 

READ (23, *) Y (5) 

READ (23, *) Y (6) 

READ (23,*) PZERO 
READ (23, 622) SSS 
READ (23, 622) SSS1 
READ (23, 622) SSS2 
READ (23, *) CONV 
READ (23, * ) MONOL 
READ (23, *) ENGACT 
READ (23, *) AREF 
READ (23, * ) TZERO 
READ (23, *)MVFLZ 
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READ (23,*) FACTOR 
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401 

READ (23, *) DX 


402 

READ (23 , * ) DENCAT 


403 

READ (23,*) EMI CRO 


404 

READ (23,*) EMACRO 


405 

READ (23, *) RMICRO 


406 

READ (23, *) RMACRO 


407 

CLOSE (UNIT = 23) 


408 

C If no external data file exists, then 

assign default values 

409 

ELSE 


410 

2119 WRITE (*,*) 


411 

WRITE (*,*)’ Default values being assigned...' 

412 

11 MFDH=100 . 0 


413 

MFDV=100 . 0 


414 

MCST=0 . 5 


415 

MCW=4 . 0 


416 

MCH=4 . 0 


417 

MCWT=0 . 25 


418 

Y (1 ) =0 . 2 


419 

Y(2)=0. 0000001 


420 

Y (3) =0 . 37 


421 

Y (4) =0 . 02 


422 

Y (5) =0.01 


423 

Y (6) =0 . 40 


424 

PZERO=l . 0 


425 

SSS='A' 


426 

SSS1='P' 


427 

SSS2='0' 


428 

CONV=2 . 5 


429 

MONOL=5 . 0 


430 

ENGACT=39700.0 


431 

AREF=123 . 4 


432 

TZERO=300 . 0 


433 

MVFLZ=0 . 25 


434 

FACTOR=50 


435 

DX=0 . 1 


436 

DENCAT=5. 0E-04 


437 

EMICRO=0 . 24 


438 

EMACRO=0.48 


439 

RMICRO=l . 2E+1 


440 

RMACRO=5 . OE+2 


441 

ENDIF 


442 

C Compute a few needed parameters: 


443 

’ Y10=Y (1) 


4 44 

Y20=Y (2) 


445 

Y30=Y (3) 


446 

Y40=Y (4) 


447 

Y50=Y (5) 


448 

Y60=Y (6) 


449 

O2ZERO=100*Y (5) 


450 

CHWP = 2 * ( MCW + MCH ) 


451 

CHSTV = CHWP / ( MCW * MCH ) 


452 

CHHD = 4 * ( MCW * MCH ) / CHWP 


453 

NOCHH = (MFDH / (MCH+2*MCST) ) 


454 

NOCHV = (MFDV / (MCW+2*MCST) ) 


455 

NOCHM = NOCHV * NOCHH 


456 

CDI=DX*FACTOR 
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9001 CONTINUE 

C Set FIRST equal to true the first time through the loop 
C so that an intial parameters may be viewed. 

FIRST = .TRUE. 

C Print Out Main Menu 
CALL LINES (17) 

WRITE (*,46) 

46 FORMAT (1H ,T20, 'LASCAT Main Menu') 


465 

MAIN = .TRUE. 



466 9011 

WRITE (*,*) 



467 

WRITE (*,*) ' 

1) 

Read in new operating parameters' 

468 

WRITE (*, *) • 

2) 

Show current operating parameters' 

469 

WRITE (*,*) • 

3) 

Change operating parameters' 

470 

WRITE (*,*) • 

4) 

Run program' 

471 

WRITE (*,*) » 

5) 

Exit program' 

472 

WRITE (*, *) 



473 

WRITE (*,*) • 


Type in number corresponding to ch 

474 

CALL LINES (12) 



475 

ASK=0 



476 

Read(*,*) ASK 



477 

IF (ASK.EQ.l) 

GO 

TO 9000 

478 

IF (ASK.EQ.2) 

GO 

TO 9002 

479 

IF (ASK. EQ. 3) 

GO 

TO 100 

480 

IF (ASK.EQ.4) 

GO 

TO 348 

481 

IF (ASK.EQ.5) 

GO 

TO 99999 


WRITE (*,*) BEEP, BEEP, BEEP 
WRITE (*, *) 

WRITE (*,*)' A number, 1 to 5, is required.' 

WRITE (*, *) 

WRITE (*,*) 

GOTO 9011 

C Output parameters to screen for review 

9002 CALL OUTPUT (0, MCST, MFDH, MFDV, MCW, NOCHV, NOCHH, MCWT, 

+ Y , MVFLZ , TZERO , PZERO , DENCAT , AREF , ENGACT , 

+ EMICRC, EMACRO, RMICRO, RMACRO, SSS, SSS1, SSS2, 

+ CONV , MONOL, DX, CDI ) 

IF (MAIN) THEN 
GOTO 9001 

ELSE 

GOTO 100 

ENDIF 

C^*»-*^*p £r£met:er Modification SBCTion**^**********'*’* *■**■*'*'*''*’'*'*■'*''*' ^^ ■*■*-:*■ •*-*-* 


4 99 

M 

O 

o 

MAIN = . FA1 

EE . 



C nr : 
^UL' 


CALL LINES (9) 



501 


WRITE (*,73) 




502 

73 

FORMAT (IB , 

720, 

' Paramei 

:er Modification Menu') 

503 


WRITE (*,*) 




504 

334 

WRITE (*,*) 




505 


WRITE (*, *) 

' 1) 

Change 

monolith physical dimensions' 

506 


WRITE (*, *) 

' 2) 

Change 

inlet gas composition' 

507 


WRITE (*, *) 

' 3) 

Change 

inlet gas volumetric flowrate' 

508 


WRITE (*,*) 

' 4) 

Change 

inlet gas temperature' 

509 


WRITE (*,*) 

' 5) 

Change 

inlet gas pressure ' 

510 


WRITE (*,*) 

’ 6) 

Change 

catalyst activation energy' 

511 


WRITE (*,*) 

' 7) 

Change 

catalyst reaction rate constant' 

512 


WRITE (*,*) 

' 8) 

Change 

catalyst active layer aensiry' 

513 


WRITE (*,*) 

' 9) 

Change macro/micro pore radius & void-fraction 
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514 WRITE (*,*) '10) Change thermal operation (Adiabatic/Isothermal) ' 

515 WRITE (*,*) 'll) Change output profile (Full/ Summary) ' 

516 WRITE (*,*) '12) Change termination variable (Conversion/Length)' 

517 WRITE (*,*) ’13) Change computation display interval' 

518 WRITE (*,*) 

519 WRITE (*,*) '14) Show current operating parameters' 

520 WRITE (*,*) '15) Save current operating parameters' 

521 WRITE (*,*) *16) Run program' 

522 WRITE (*,*) ’17) Return to main menu' 

523 WRITE (*,*) ’18) Exit program’ 

524 WRITE (*,*) 

525 WRITE (*,*) ' Type in number corresponding to choice above' 

526 WRITE (*,*) 

527 ASK=0 

528 READ (*, *) ASK 

529 IF (ASK.EQ. 1) GOTO 101 

530 IF (ASK.EQ.2) GOTO 110 

531 IF (ASK.EQ. 3) GOTO 120 

532 IF (ASK.EQ.4) GOTO 130 

533 IF (ASK.EQ.5) GOTO 140 

534 IF (ASK.EQ.6) GOTO 150 

535 IF (ASK. EQ. 7) GOTO 160 

536 IF (ASK.EQ.8) GOTO 170 

537 IF (ASK.EQ.9) GOTO 180 

538 IF (ASK.EQ. 10) GOTO 200 

539 IF (ASK.EQ. 11) GOTO 210 

540 IF (ASK.EQ. 12) GOTO 220 

541 IF (ASK.EQ. 13) GOTO 230 

542 IF (ASK.EQ. 14) GOTO 9002 

543 IF (ASK.EQ. 15) GOTO 11000 

544 IF (ASK.EQ. 16) GOTO 348 

545 IF (ASK.EQ. 17) GOTO 9001 

546 IF (ASK.EQ. 18) GOTO 99999 

547 WRITE (*,*) BEEP, BEEP, BEEP 

548 WRITE (*,*) 

549 WRITE (*,*)' A number, 1 to 18, is required.' 

550 WRITE (*,*) 

551 GOTO 334 

552 c************* PARAMETER INPUT SECTION************************ 

553 C INPUT MONOLITH PHYSICAL DIMENSIONS 

554 101 WRITE (*,*) 

555 MCST=MCST*2 . 0 

556 WRITE (*,*)' Current monolith horizontal ' 

557 WRITE (*,1075) MFDE 

558 WRITE (*,*) 

559 1075 FORMAT (' face dimension (mm) = ’ , F8 . 3 ) 

560 1136 WRITE (*,*)' Input new monolith horizontal face dimension.' 

561 WRITE (*,*) ' Value must be >=10mm and <=1000mm. ' 

562 READ (*, *) MFDH 

563 IF ( (MFDH. LT. 10.0) .OR. (MFDH. GT. 1000.0) ) GOTO 1136 

564 WRITE (*, *) 

565 WRITE (*,*) 

566 WRITE (*,*)' Current monolith vertical' 

567 WRITE (*,1076) MFDV 

558 1076 FORMAT ( ' face dimension (mm) = ’ , F8 . 3 ) 

569 WRITE (*,*) 

570 WRITE (*,*) 
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571 

572 

573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 

615 

616 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 
627 


1137 WRITE (*,*)’ Input new monolith vertical face dimension.' 

WRITE (*,*)' Value must be >=10mm and <=1000mm. ' 

READ (*, * ) MFDV 

IF ( (MFDV. LT. 10.0) .OR. (MFDV. GT. 1000.0) ) GOTO 1137 

102 WRITE (*, *) 

WRITE (*, *) 

WRITE (*,*)' Current monolith channel total' 

WRITE (*,1077) MCST 

1077 FORMAT (’ wall thicknes (mm) = ',F8.3) 

WRITE (*, *) 

WRITE(*, *) 

WRITE (*,*) ’ Input the monolith channel total wall' 

WRITE (*,*)' thickness (mm) . (Including active layer)' 

READ (*, *) MCST 

IF (MCST. LE. 0.0) GO TO 102 

IF ( (MCST . GE . MFDH ) . OR . (MCST . GE . MFDV) ) THEN 

WRITE (*,*)' Wall thickness must be less than the monolith' 
WRITE (*,*)' face dimension! Right? Let us start over.' 

GOTO 101 
ENDIF 

MCST = MCST / 2.0 

103 WRITE (*,*) 

WRITE (*,*)' Current square channel' 

WRITE (*,1078) MCW 

1078 FORMAT (' opening dimension (mm) = ',F8.3) 

WRITE (*, * ) 

WRITE (*,*)' Input new monolith square' 

WRITE (*,*)' channel opening dimension (mm) . ' 

READ (*, *) MCW 

IF (MCW. LE. 0.0) GO TO 103 

IF ( (MCW. GE. MFDH) .OR. (MCW. GE. MFDV) ) THEN 

WRITE (*,*) ' Channel opening must be less than the monolith' 
WRITE (*,*)' face dimension! Right? Let us start over.' 

GOTO 101 
El© IF 
MCH = MCW 

104 WRITE (* , * ) 

WRITE (*,*)' Current active layer thickness' 

WRITE (*,1079) MCWT 

1079 FORMAT (' thickness (mm) * ’,F8.3) 

WRITE (*,*) 

WRITE (*,*) 

WRITE (*,*)' Input monolith channel active layer ' 
write (*,*)' thickness (mm). Active layer cannot' 

WRITE (*,1080) MCST* 2. 00 

1080 FORMAT (' be greater than ',F8.3, ' (mm),') 

WRITE (*,*) ' (1/2 total wall thickness) ' 

READ (*, *) MCWT 

IF ( (MCWT. LE.0) .OR. (MCWT. GT. MCST) ) GO TO 104 
WRITE (*,*) 

WRITE (*,*) 

C Perform monolith parameter calculations 
8900 CHWP = 2 * ( MCW + MCH ) 

CHSTV = CHWP / ( MCW * MCH ) 

NOCHH = (MFDH / (MCH+2*MCST) ) 

NOCHV = (MFDV / (MCW+2*MCST) ) 

CHHD =4.0* (MCW*MCH) /CHWP 



I 
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628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 


NOCHM = NOCHV * NOCHH 
GO TO 100 

Input Gas composition 


C 

C 

C 

110 


668 

669 

670 

671 

672 

673 

674 

675 

676 

677 
676 

679 

680 
681 
682 

683 

684 


Define Mole fractions of species Y(i) 

The index of Y(i) correspond to the species. 
Species : l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2 . 


Write (*, *) 
WRITE(*,*) 
Write (*, *) 
Write!*,*) 
Write!*,*) 
WRITE!*,*) 


Input section for gas composition:' 


WRITE!*,*) 


3231 

3232 


The mixture can consist of only the following' 
species: He, Ar, C02, CO, 02, N2. All mole ' 
fractions must be greater than zero. The sum' 
of the mole fractions must be =1' 

Current mole fractions: ' 

WRITE!*, 3231) Y(l) , Y (2) , Y (3) 

WRITE!*, 3232) Y (4) , Y (5) , Y (6) 

FORMAT!' He = ',F6.4,' Ar = ',F6.4,' 

FORMAT!' CO = ' , F6 . 4, ' 02 = \F6.4,' 

WRITE!*,*) 

Write!*,*) 'Input gas mole fraction for He' 

Read!*,*) Y(l) 

If (Y(l) .EQ.0.0) Y (1) =lE-7 


C02 = ' ,F6. 4) 
N2 = ' ,F6.4) 


649 

Write (*, *) 

'Input gas mole fraction for Ar' 


650 

Read!*, * 

) Y (2) 


651 

If (Y (2) 

• EQ. 0.0) Y (2) =lE-7 


652 

Write (*, *) 

'Input gas mole fraction for C02' 


653 

Read!*, * 

) Y (3) 


654 

If (Y (3) 

•EQ.0.0) Y (3) =lE-7 


655 

WRITE!*,*) 



656 

Write!*,*) 

'Input gas mole fraction for 02' 


657 

WRITE (*, *) 



658 

Write (*, *) 

'Gas mole fraction for CO will be 

set to twice that 

659 

Write (*, *) 

'of 02. This is because CO and 02 

are presumed to' 

660 

Write (*, *) 

'be produced from the dissociation 

of C02 . ' 

661 

Write (*, *) 

'Therefore, the stoichiometric 2:1 

ratio applies . ' 

662 

Write!*,*) 

'Also, the mole fraction of 02 is 

limited to a ' 

663 

Write (*, * ) 

'maximum of 0.03 (3.0%). Typically a C02 lasers ' 

664 

write (*, *) 

'output is degraded at a 02 mole fraction of 0.01’ 

665 

Write (*, *) 

' (1.0%) . ' 


666 

Read!*, * 

) Y (5) 


667 

IF (Y (5) 

.LE.0.0) THEN 



Write (*, *) 
Write!*,*) 

GOTO 110 
END IF 

IF (Y (5) .LT.1E-7) THEN 


There must be 
Put in a mole 


02 to convert 
fraction > O' 


!!', BEEP, BEEP 


no problem! Put in 
BEEP, BEEP 


WRITE!*,*) 

Write!*,*) 

Write (*, *) 

GOTO 110 
END IF 

IF (Y(5) .GT. 0.0300001) THEN 
Write!*,*) ' Awful high 
Write!*,*) 

Write (*,*) 

Write (*, *) 

GOTO 110 
ENDIF 


Awful low 02 concentration. The problem is' 


a mole fraction >= IE-7' 


02 concentration. The laser' 
may explode! Shut it off immediately' 
BEEP, B EEP 

Put in a mole fraction <=0.03 ' 
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685 


Y (4)=2.0*Y (5) 


686 


02ZERO=Y (5) *100 


687 


Y (6) = 1 . 0-Y (5) -Y (4) -Y (3) -Y (2) -Y (1) 


688 


IF ( (Y (6) .LT.0.0) .OR. (Y(6) .GT.1.0) ) GOTO 110 


689 


IF ( (Y (6) .LT.1E-7) .AND. (Y (6) .GE.O) ) Y(6)=lE-7 


690 


WRITE (*,*) 


691 

C Check if restrictions are met, If not, force user to 

reinput values 

692 


Do 112 1=1,5 


693 


IF ( (Y (I) .LT.0) .or. (Y (I) .GT.1.0) ) THEN 


694 


WRITE (*,*) BEEP, BEEP 


695 


WRITE (*,*)' Invalid mole fractions. Try again.' 

696 


WRITE (*,*)' Restrictions are not met. 

Read below' 

697 


Write (*,*) 


698 


GOTO 110 


699 


ENDIF 


700 

112 

Continue 


701 


WRITE (*,2323) Y(6) 


702 

2323 

FORMAT (' The gas mole fraction for N2 = ',F6.4) 


703 


Y10=Y (1) 


704 


Y20=Y (2) 


705 


Y30=Y (3) 


706 


Y 4 0=Y ( 4 ) 


707 


Y50=Y (5) 


708 


Y60=Y (6) 


709 


WRITE (*, *) 


710 


WRITE (*,*) 


711 


WRITE (*,*) * — > Hit Return when Finished Viewing <=* 

712 


WRITE (* , * ) 


713 


READ (*, 621 ) S 


714 


WRITE (*, * ) 


715 


WRITE (*, * ) 


716 


GO TO 100 


717 

C Input Inlet gas flowrate 


718 

120 

WRITE (*,*) 


719 


WRITE (* , * ) 1 Current inlet volumetric* 


720 


WRITE (*, 1081 ) MVFLZ 


721 

1081 

FORMAT ( 1 flowrate = , / F8.3 / l (1/s) . 1 ) 


722 


WRITE(*,*) 


723 


WRITE (*,*) * Input monolith inlet volumetric flowrate.* 

724 


READ (*, *) MVFLZ 


725 


IF (MVFLZ. LE. 0.0) GO TO 120 


726 


GO TO 100 


727 

C Inn- 

ut Inlet gas temperature 


728 

130 * 

WRITE (*, *) 


729 


WRITE (*,*)' Current inlet aas ' 


730 


WRITE (*,1082) TZERO 


731 

1082 

FORMAT ( ' temperature = ' , F8 . 3, 1 (K) , ' ) 


732 


WRITE (*,*) 


733 


WRITE (*,*) 'Input monolith inlet gas temperature. 

I 

734 


WRITE (*,*) 'Temperature must be >200K and <700K. ' 


735 


READ (*, *) TZERO 


736 


IF ( (TZERO. EQ.O) .OR. (TZERO.LT. 200) ) GO TO 130 


737 


IF (TZERO. GT. 700) GO TO 130 


738 


GOTO 100 


739 

C Input Inlet Gas Pressure 


740 

140 

WRITE (*,*) 


741 


PZERO=PZERO*101 . 325 
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742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

771 

772 

773 

774 

775 

776 

777 

778 

779 

780 

781 

782 

783 

784 


786 

787 

788 

789 

790 
7 91 

792 

793 

794 

795 

796 

797 

798 


WRITE (*,*)' Current inlet gas ' 

WRITE (*,1083) PZERO 

1083 FORMAT ( ' pressure = ’,F8.3,' (kPa) ' ) 

WRITE (*,*) 

WRITE (*,*) ' Inlet gas pressure must be between' 

WRITE (*, *) ' 0.1 and 350 kPa' 

WRITE (*,*)' Input new inlet gas pressure.' 

READ (*, *) PZERO 

IF ( (PZERO. LT. 0.1) .OR. (PZERO. GT. 350.0) ) GOTO 140 
C convert pressure back to atm. 

PZERO=PZERO/101 . 325 
GOTO 100 

C Input catalyst activation energy 
150 WRITE (*, * ) 

WRITE (*,1084) ENGACT 

1084 FORMAT ( ' Current activation energy = ',F11.3, ' (J/mol),') 
WRITE (*,*)' Maximum activation energy (J/mol) = 300000' 
WRITE (*,*) 

WRITE (*,*)' Catalyst activation energy ( J/mol ) ' 

READ (*, *) ENGACT 

IF ( (ENGACT. LE. 0.0) .OR. (ENGACT. GT. 300000) ) GOTO 150 
GO TO 100 

C Input catalyst reaction rate constant 
160 WRITE (*,*) 

WRITE (*,*)' Current catalyst reaction rate' 

WRITE (*,1085) AREF 

1085 FORMAT (’ constant [@ 298 K] = ',F11.3, ' (mm"3/gcat-s) , ' ) 
WRITE (*,*) 

WRITE (*,*)' Input new catalyst reaction' 

WRITE (*,*)' rate constant at 298 K (mm / '3/gcat-s) ' 

WRITE (*,*) ' Rate must be greater than zero' 

READ ( * , * ) AREF 
IF (AREF. LE. 0.0) GOTO 160 
GO TO 100 

C Input catalyst active layer density 
170 WRITE (*,*) 

write (*,*)' C urrent density of catalyst' 

WRITE (*,1086) DENCAT 

1086 FORMAT ( ' active layer = ',E10.4, ’ (g/mm~3),’) 

WRITE (* , * ) 

WRITE (*,*)' Input density of catalyst active layer (g/mm^3) ' 
WRITE (*,*) 'Density must be greater than zero' 

READ (*, *) DENCAT * 

IF (DENCAT. LE. 0.0) GOTO 170 
WRITS (*, *) 

GOTO 100 

C Input micro/macropore radii and void-fractions 
180 * WRITE (*,*) 

WRITE (*,*)' Current average micropore ' 

WRITE (*, 1087) RMICRO 

1087 FORMAT (' radius = ',F9.4,' (nm) , ’ ) 

WRITE (*,*)' Input new average micropore radius (nm) ' 

WRITE (*,*) ' Micro radius must be greater than zero' 

READ ( * , * ) RMICRO 
WRITE (*,*) 

WRITE (*,*)' Current average macropore' 

WRITE (*,1088) RMACRO 


1 



799 1088 

800 
801 
802 

803 

804 

805 

806 

807 

808 

809 

810 
811 
812 

813 190 

814 

815 1089 

816 1090 

817 

818 

819 

820 
821 
822 

823 

824 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834 

835 

836 

837 

838 

839 

840 

841 C Ira 

842 200 

843 

844 

845 

846 

847 

848 

849 

850 

851 

852 

853 

854 


FORMAT ( ' radius = ',F9.4,' (nm) , ' ) 

WRITE (*,*) ’ Input new average macropore radius (nm) ' 

WRITE (*,*) ' Macro radius must be greater than zero' 

READ (*, *) RMACRO 
WRITE (*, *) 

IF ( (RMICRO.LE. 0.0) .OR. (RMACRO. LE. 0.0) ) GO TO 180 
IF (RMICRO.GT. RMACRO) THEN 

WRITE (*,*)' Hey Dude, What's the deal?’ 

WRITE (*,*)’ You just entered macropore radius which is' 
WRITE (*,*) ' smaller than the micropore radius! ! ' 

WRITE (*,*)' Try again, big guy !!' 

WRITE (*,*) 

GOTO 180 

ENDIF 

WRITE (*,1089) EMICRO 
WRITE (*,1090) EMACRO 

FORMAT (' Current void-fraction as micropores ', F6.3) 

FORMAT (' Current void-fraction as macropores ', F6.3) 

WRITE (*,*) 

WRITE (*,*) ' The sum of micro and macropore void-fractions' 
WRITE (*,*) ' must be less than one (1) ' 

WRITE (*,*) ' Void fractions must be >=0 and <=1' 

WRITE (*,*) 

WRITE (*,*)' Input void-fraction as micropores ' 

READ (*, *) EMICRO 

WRITE (*,*)' Input void-fraction as macropores ' 

READ (*, *) EMACRO 
WRITE (*, *) 

IF ( (EMICRO.LT. 0.0) .OR. (EMACRO. LT. 0 . 0) ) THEN 

WRITE (*,*)' You can't have a negative void fraction’ 
WRITE (*,*) 

WRITE (*, *) 

GOTO 190 

ENDIF 

IF ( EMI CRO+ EMACRO . GT . 1 . 0 ) THEN 

WRITE (*,*) 1 The sum of micro and macropcre void' 

WRITE (*,*) ' fractions must be less than one (1) ’ 

WRITE (*,*) 

GOTO 190 

ENDIF 
WRITE (*,*) 

GOTO 100 

>ut thermal opearation: adiabatic cr isothermal 
IF (SSS.EQ.*'!' ) THEN 

WRITE (*,*)' Current operation is isothermal' 

ELSE 

WRITE (*,*)' Current operation is adiabatic' 

END IF 
WRITE (*,*) 

WRITE (*,*) 'Isothermal (I) or Adiabatic (A) operation?' 

WRITE (*,*) 'Type the letter I or the letter A. HIT RETURN' 

READ (*,621) SSS 
IF (SSS.EQ. 'i') SSS = 'I' 

IF (SSS.EQ. 'a' ) SSS = 'A' 

IF ( (SSS.NE. 'I' ) .AND. ( SSS.NE.'A')) GOTO 200 
GO TO 100 


855 C Input output profile type: full or summary 
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856 

857 

858 

859 

860 
861 
862 

863 

864 

865 

866 

867 

868 

869 

870 

871 

872 

873 

874 

875 

876 

877 

878 

879 

880 
881 
882 

883 

884 

885 

886 

887 

888 

889 

890 
831 

892 

893 

894 

895 

896 
637 

898 

899 

G 

“uO 
Qni 
u — 

9C2 

903 

904 

905 

906 

907 

908 

909 

910 

911 

912 


210 WRITE (*,*) 

WRITE (*,*)' Full concentration and temperature profile (P) or' 
2 ' summary (S) ? ' 

WRITE (*,*) 'Type the letter P or the letter S. HIT RETURN. ' 
READ (*,621) SSS1 
IF (SSS1.EQ. 'p') SSS1 = ’P’ 

IF (SSS1.EQ. ’s') SSS1 = 'S' 

IF ((SSS1.NE. 'P') .AND. ( SSS1.NE. ’S’) ) GO TO 210 
GO TO 100 

C Input computation loop completion criteria 
220 WRITE (*,*) 

WRITE (*,*) ' Integration can be stopped by either specifying' 
WRITE (*,*)' an amount of 02 conversion (O) , or by specifying' 
WRITE (*,*)' a monolith length (L) ' 

WRITE (*,*) 

IF (SSS2 .EQ. ’O') THEN 

WRITE (*,*)' Currently, Integration will be stopped' 

WRITE (*,*)' by specifying an amount of 02 conversion (O).' 
ELSE 

write (*,*)' Currently, Integration will be stopped' 

WRITE (*,*) ' by specifying a monolith length (L) . ' 

END IF 
WRITE (*,*) 

WRITE (*,*) ' Type the letter 0 or the letter L. HIT RETURN. ' 
READ (*,621) SSS2 
IF(SSS2.EQ. ’o') SSS2 = 'O' 

IF (SSS2 . EQ . ' 1 ' ) SSS2 = 'L' 

IF ( (SSS2.NE. 'O' ) .AND. { SSS2.NE. ' L ' ) ) GO TO 220 
IF (SSS2.EQ. 'O' ) THEN 
225 IF (MONOL.EQ.O) THEN 

WRITE (*,*) 

WRITE (*,1091) CONV 

1091 FORMAT (' Current % oxygen conversion « ',F8.3) 

END IF 
WRITE (*, *) 

WRITE (*,*)' Input desired % oxygen conversion' 

WRITE (*, * ) ' between 0 and 100%.' 

READ (*, *) CONV 

IF ( (CONV. GT. 100.0) .OR. (CONV. LT.0) ) GOTO 225 

ENDIF 

227 IF (SSS2.EQ. 'L' ) THEN 
WRITE(’,*) 

WRITE (*,1092) MONOL 

1092 FORMAT (' Current desired mcnolith length = ' ,F8.3, 1 (mm)') 

WRITE (*, *) 

WRITE (*,*) ’ Input desired monolith lencrth (5-4000mm) ' 
WRITE (*, *) 

READ (*, *) MONOL 
IF (MONOL. LT. 5.0) GOTO 227 
IF (MONOL. GT. 4000.0) GOTO 227 

ENDIF 
GO TO 100 

C Input computation loop (intecration) steD size *** 

230 WRITE (*, *) 

WRITE (*,1093) CD I 

1093 FORMAT (' Current computation display interval ', F6.3, ' (mm)') 

WRITE (*,*)' Minimum computation display interval = 0.001 mm.’ 
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913 

914 

915 

916 

917 

918 

919 

920 

921 

922 

923 

924 

925 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

940 

941 

942 

943 

944 

945 

946 

947 

948 

949 


WRITE (*,*)' Maximum computation display interval = 100 mm.' 
WRITE (*,*)' Values > 1 will be rounded to the nearest integer' 
WRITE (*,*) 

WRITE (*,*)' Input desired computation display interval (mm).' 
WRITE (*, *) 

WRITE (*,*) 

READ (*, *) CD I 
IF (CDI.GT.100) GOTO 230 
IF (CDI.GT.1.0) CDI=INT (CDI) 

IF (CDI. LT. 0.0009999) GOTO 230 
WRITE (*,*) 

WRITE (*,*) 

GO TO 100 

C END of parameter input section 
C Save present operationg parameters to file 
11000 WRITE (*,*) 

C Give user option to save current operating parameters. 

WRITE (*,*) 

WRITE (*,*) ' Would you like to save the current operating' 
WRITE(*,*)' parameters to a New(N) or Existing(E) file?' 

WRITE (*,*) 

READ (*,621) CHI 

IF( (CH1.NE. 'E' ) .AND. (CH1.NE. 'e' ) .AND. (CH1.NE. 'N') 

* .AND. (CH1.NE. 'n' ) ) GOTO 11000 

IF ((CH1.EQ. 'E') .OR. (CH1.EQ. ’ e ' ) ) THEN 
WRITE (*, *) 

WRITE (*,*) 'What's the name of the EXISTING’ 

WRITE (*,*) 'parameter file to be opened?' 

READ (*,621) FILNAM 

OPEN ( UNIT=24, FILE=FILNAM, STATUS= ' OLD ' ) 

ELSE 

WRITE (*,*) 

WRITE (*,*) 'What' s the name of the NEW' 

WRITE (*,*) 'parameter file to be opened?' 

READ (*,621) FILNAM 

OPEN ( UNIT— 2 4 , FILE=FILNAM, STATUS= ' NEW ' ) 

ENDIF 


950 C Write default parameters to parameter file. 


951 

WRITE (24, *)MFDH 

952 

WRITE (24, *)MFDV 

953 

WRITE (24,*) MOST 

954 

WRITE (24, *)MCW 

0 ^ K 

^ W nj 

WRITE (24,*) MCE 

956 

WRITE (24, *)MCWT 

957 

WRITS (24, *)Y(1) 

958 

WRITE (24, *) Y (2) 

959 

WRITE (24, *) Y (3) 

960 

WRITE (24, *) Y (4) 

961 

WRITE (24, *) Y (5) 

962 

WRITE (24, *) Y (6) 

963 

WRITE (24,*) PZERO 

964 

WRITE (24, 622) SSS 

965 

WRITE (24, 622) SSS1 

966 

WRITE (24, 622) SSS2 

957 

WRITE (24, *) CONV 

968 

WRITE (24, *)MONOL 

969 

WRITE (24, *)ENGACT 
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970 

971 

972 

973 

974 

975 

976 

977 

978 

979 

980 

981 

982 

983 

984 

985 

986 

987 

988 

989 

990 

991 

992 

993 

994 

995 

996 

997 

998 

999 
1000 
1001 
1002 

1003 

1004 

1005 

1006 

1007 

1008 

1009 

1010 
*1 m i 

-k.- 

1012 

1013 

1014 

1015 

1016 

1017 

1018 

1019 

1020 
1021 
1022 

1023 

1024 

1025 

1026 


WRITE (24, *)AREF 
WRITE (24, * ) TZERO 
WRITE (24, *)MVFLZ 
WRITE (24,*) FACTOR 
WRITE (24, *)DX 
WRITE (24,*) DENCAT 
WRITE (24, *)EMICRO 
WRITE (24, *)EMACRO 
WRITE (24, *)RMICRO 
WRITE ( 2 4 , * ) RMACRO 
CLOSE ( UNIT = 24) 

GOTO 100 

C COMPUTE GAS AVG. MOLECULAR WEIGHT AND DENSITY * 

348 AMW = FXN4 (Y) 

R=8 . 3144 

DENGAS=AMW* (PZERO/9896.0) /R/TZERO 
C Parameters needed for thermal cond. calcuation 
C Read in viscosity (VISC) values into matrix. Units for 
C VISC are centipoises. The values are read off a nomagraph 

C in Perry's Chemical Engineering Handbook, 5th ed. , p3-210. 

C The first index of VISC corresponds to the species. 

C Species : l=He , 2=Ar , 3=C02 , 4=CO, 5=02 , 6=N2 . 

C The second index of VISC corresponds to the temperature. 

C Temperature : 2=200K, 3=300K, 4=400K, 5=500K, 6=600K, 7=700K 

DATA VISC (1,2) , VISC (1, 3) , VISC (1,4) ,VISC (1, 5) , VISC (1, 6) , VISC (1,7) 
+/.014, .0188, .022, .027, .031, .034/ 

DATA VISC (2,2) , VISC (2, 3) , VISC (2, 4) ,VISC (2, 5) , VISC (2, 6) ,VISC (2,7) 
+/.0164, .0223, .0278, .0327, .0376, .041/ 

DATA VISC (3, 2) ,VISC(3,3) ,VISC(3,4) ,VISC(3,5) , VISC (3, 6) ,VISC(3,7) 
+/.0105, .0147, .019, .023, .027, .031/ 

DATA VISC (4, 2), VISC (4, 3), VISC (4, 4), VISC (4, 5), VISC (4, 6), VISC (4, 7) 
+/. 0138, .018, .022, .0255, .029,-0325/ 

DATA VISC (5, 2) ,VISC(5,3) , VISC (5, 4) , VISC (5, 5) ,VISC(5,6) ,VISC(5,7) 
+/.0155, .02, .0247, .0286, .0327, .0368/ 

DATA VISC (6, 2) ,VISC(6,3) ,VISC(6,4) ,VISC(6,5) ,VISC(6,6) ,VISC(6,7) 
+/.013, .0175, .0218, .026, .029, .0332/ 

MVFLO=MVFLZ 

C CONVERT MVFLO FROM liters/s TO mm^3/s 
MVFLO = MVFLO *1.0£+06 
C COMPUTE MOLAR FLOWRATE OF GAS. (M/s) 

MOLFLO = MVFLO*DENGAS * AMW 
C COMPUTE CHANNEL FLOW RATE AND CHANNEL VELOCITY 
CHFLO = MVFLO / NOCHM 
CHGVEL = CHFLO / ( MOW * MCH) 

349 CONTINUE 

C INTIALIZE PHYSICAL PROPERTIES PARAMETERS 
1041 HEATRX=5. 64E-r05 
TRZF=298 . 0 
PRESS = PZERO 
TGX= TZERO 
TSX= TZERO 

CGX= 4 . 092E-08 * (O2ZERO/100 . 0) * (298. 0/TGX) * (PRESS/1 . 0) 

FGX=CGX* CHFLO 
FZERO=FGX 

FGXEND=F ZERO* (1.0- CONV/100 . 0 ) 

S = CHAR (9) 
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1027 

1028 

1029 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 

1038 

1039 

1040 

1041 

1042 

1043 

1044 

1045 

1046 

1047 

1048 

1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 

1058 

1059 

1060 
1061 
1062 

1063 

1064 

1065 

1066 
1067 
1066 
-i w'C-/ 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1078 

1079 

1080 
1081 
1082 
1083 
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SS = CHAR (13) 

X=0 . 0 

CMCWT= ( (MCH+2*MCWT) * (MCW+2*MCWT) - (MCH*MCW) ) /2/ (MCH+MCW) 

CONSTA=0 . 0 

DCGX=0 . 0 

DFGX=0 . 0 

ICNT = 0 

C Initialize variables for Runge Kutta Integration method 
XNEW(l) = FGX 
XNEW (2) = PZERO 
XNEW (3) = TZERO 

C************* START OF STEP SIZE CHECK LOOP ******************** 
1060 DO 8001 P=0, 2 
DX=lE-6 

DENGAS=AMW* (PRESS/9896.0) /R/TGX 

C COMPUTE MONOLITH FLOW RATE, CHANNEL FLOW RATE AND CHANNEL VELOCITY 
MVFLO = MOLFLO / (AMW * DENGAS) 

CHFLO = MVFLO / NOCHM 
CHGVEL = CHFLO / ( MCW * MCH) 

C COMPUTE BULK DIFFUSIVITY 
DAB = FXN3 (Y, TGX) 

C COMPUTE EFFECTIVE DIFFUSIVITY 

DABEFF= FXN1 (DAB, Y (5) , PRESS, RMICRO, RMACRO, 

+ TSX, EMACRO, EMICRO) 

C COMPUTE THERMAL CONDUCTIVITY 
T=TGX 

C COMPUTE VISCOSITY OF COMPONENT GASES AT TGX 
CALL VIS (TGX, VISC, VISCT) 

C COMPUTE THERMAL CONDUCTIVITY OF COMPONENT GASES AT TGX 
CALL TCON (TGX, TCONT) 

C COMPUTE VISCOSITY AND THERMAL CONDUCTIVITY OF MIXTURE 
CALL PROP (Y, TCONT, VISCT, VISCTC,K) 

C COMPUTE HEAT CAPACITY OF GAS MIXTURE AT TGX 
CPGAS=FXN2 (Y, TGX) 

KC = 2. 976 * DAB / CHHD 
H = 2. 576 * K / CHHD 

A = AREF * DEXP ( -ENGACT / R * ( 1 / TSX -1 / TREF) ) 

THIELE = CMCWT * ( A * DENCAT / DABEFF) ** 0.5 

N = DTANH ( THIELE ) / THIELE 

C0NST1= N*A*CMCWT*DENCAT 

CSX = CGX * ( KC) / (KC+CONST1) 

TSX =TGX + KEATRX*CSX*CONSTl/E 

C FIRST TIME THROUGH COMPUTE 'A' TO USE FOR TSX, CSX CALCULATION 
IF (FIRST. EQ. .TRUE. ) THEN 

A = AREF * DEXP ( -ENGACT / R * ( 1 / TSX -1 / TREF) ) 

THIELE = CMCWT * ( A * DENCAT / DABEFF) **0.5 

N = DTANH ( THIELE ) / THIELE 

CONSTl= N*A*CMCWT*DENCAT 

CSX = CGX * ( KC) / (KC+CONST1 ) 

TSX =TGX + KEATRX* CSX* CONST1 / H 
IF (SSS.EQ. ’ I * ) TSX=TGX 
IF (P.EQ.O) THEN 
XZ=0 

CONVXZ=0 . 0 

CGXZ=CGX 

CSXZ=CSX 
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1084 

1085 

1086 

1087 

1088 

1089 

1090 

1091 

1092 

1093 

1094 

1095 

1096 

1097 

1098 

1099 
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1101 
1102 

1103 

1104 

1105 

1106 

1107 

1108 

1109 

1110 
1111 
1112 

1113 

1114 

1115 

1116 

1117 

1118 

1119 

1120 
1121 
1122 

1123 

1124 

1125 

1126 

1127 

1128 

1129 

1130 
1 1 v\ 

1132 

1133 

1134 

1135 

1136 

1137 


TSXZ=TSX 

ENDIF 

ENDIF 

RE=1000 . 0*CHHD*CHGVEL*DENGAS /VISCTC 
CONVX = 100.0* (CGXZ-CGX*TGX/TZERO*PZERO/PRESS) /CGXZ 
C Reset variables each time Runge Kutta Subroutine is called 
7889 X0(1) = XNEW(l) 

XO (2) = XNEW ( 2 ) 

XO (3) = XNEW (3) 

C Call Runge Kutta Subroutine: 

CALL RK4 (DX,X,XO,XNEW) 

C PARMETER UPDATE 

FGX = XNEW (1 ) 

PRESS = XNEW (2) 

IF (SSS.EQ.’I’) THEN 
TGX = TZERO 

ELSE 

TGX = XNEW (3) 

ENDIF 

CGX = FGX/CHFLO 
X = X + DX 

Y (5) =Y50* (100 . O-CONVX) /100 . 0 
Y (4)=Y40-2* (Y50-Y (5) ) 

Y(3)=Y30+2* (Y50-Y (5) ) 

Y (1 ) =Y10 
Y (2) =Y20 
Y (6) =Y60 

CONST7 = Y (1) + Y (2) + Y(3) + Y(4)+ Y(5)+ Y(6) 

DO 33446 JJ=1 , 6 

Y(JJ)=Y(JJ) /CONST7 
33446 CONTINUE 
8001 CONTINUE 

0************* END OF STEP SIZE CHECK LOOP ******************** 


C CALCULATE STEP SIZE AND DISPLAY FACTOR USING COMPUTATIONAL 
C DISPLAY FACTOR AND CONVX CALCULATED IN STEP SIZE CHECK LOOP 
DX=lE-6*0 . 04 /CONVX 
IF (DX.GT.5) DX=5 
IF (DX.LT. 0 . IE-4) DX=0.1E-4 

IF (CDI.LT.DX) DX=CDI 
2345 FACTOR=INT ( CEI /DX ) 

IF (FACTOR. GE. 250) THEN 
CDI=CDI*250/FACTOR 
FACTOR=250 

ENDIF 

DX=CDI /FACTOR 

C INITIALIZE PARAMETERS FOR MAIN COMPUTATIONAL LOOP 
PRESS = PZERO 
TGX= TZERO 
TSX= TZERO 

CGX= 4 . 092E-08 * (O2ZERO/100 . 0) * (298 . 0/TGX) * (PRESS/1 . 0) 
FGX=CGX*CHFLO 


1138 FZERO=FGX 

1139 FGXEND=FZERO* (1.0 - CONV/100.0) 

1140 CGXTMP=CGX 
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1141 

1142 

1143 

1144 

1145 

1146 

1147 

1148 

1149 

1150 

1151 

1152 

1153 

1154 

1155 

1156 

1157 

1158 

1159 

1160 
1161 
1162 

1163 

1164 

1165 

1166 

1167 

1168 

1169 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 

1180 
1181 
1182 


1184 

1185 

1186 

1187 

1188 

1189 

1190 

1191 

1192 

1193 

1194 

1195 

1196 

1197 


X=0 . 0 

CONSTA=0 . 0 
DCGX=0 . 0 
DFGX=0 . 0 
ICNT = 0 
TGNEW=TZERO 
MVFLOMVFLZ 
Y (1) =Y10 
Y(2)=Y20 
Y(3)=Y30 
Y ( 4 ) =Y 4 0 
Y (5) =Y50 
Y ( 6 ) = Y 6 0 

C Initialize variables for Runge Kutta Integration method 
XNEW(l) = FGX 
XNEW (2) = PZERO 
XNEW (3) = TZERO 
JUMP=0 
ONCE=0 


0************* START OF MAIN COMPUTATIONAL LOOP ******************** 
1059 DO 8000 P=0, 100000 

IF ( (DCGX.GT.0.0) .OR. (CGX.LT.0.0) ) THEN 

WRITE (*,*)’ Something is wrong, conversion is proceeding' 
in the wrong direction or the concentration' 
of 02 has dropped below zero. ' 

Try a lower computation display increment. ' 

Try a lower gas inlet temperature or flowrate. ' 


WRITE (*,*) 
WRITE (*, *) 
WRITE (*,*) 
WRITE (*,*) 
WRITE (*, *) 
WRITE (*,*) 
WRITE (*,*) 
WRITE (*,*) 
GOTO 17000 
ENDIF 


Change in cone of 02 per step(DCGX) = ',DCGX 
02 concentration (CGX) = ',CGX 


C Check to see if conversion has been reached. 

IF (SSS2.EQ. 'O' ) THEN 

IF ( (CONVX.GE.CONV) ) GOTO 9900 

ELSE 

IF (X.GT. MONOL) GOTO 9900 
ENDIF 

C COMPUTE. GAS AVG. MOLECULAR WEIGHT AND DENSITY 
6250 AMW = FXN4 (Y) 

DENGAS=AMW* (PRESS/9896.0) /R/TGX 

C COMPUTE MONOLITH FLOW RATE, CHANNEL FLOW RATE AND CHANNEL VELOCITY 
MVFLC = MCLFLO / (AMW * DENGAS) 

CKFLO = M7FL0 / NOCHM 
CHGVEL = CHFLO / ( MCW * MCH) 

C COMPUTE BULK DIFFUSIVITY OF 02 
DAB = FXN3 (Y, TGX) 

C COMPUTE EFFECTIVE DIFFUSIVITY OF MIXTURE AT TSX 
DAEEFF= FXN1 (DAB, Y (5) , PRESS, RMICRO, RMACRO, 

+ TSX, EMACRO, EMICRO) 

IF (MOD ( CONSTA, 4) .GT. 0.0 ) GO TO 6669 


T=TGX 

C COMPUTE VISCOSITY OF COMPONENT GASES AT TGX 
CALL VIS (TGX, VISC, VISCT) 
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1198 C COMPUTE THERMAL CONDUCTIVITY OF COMPONENT GASES AT TGX 

1199 CALL TCON ( TGX , TCONT ) 

1200 C COMPUTE VISCOSITY AND THERMAL CONDUCTIVITY OF MIXTURE 

1201 CALL PROP (Y, TCONT, VISCT, VISCTC, K) 

1202 6666 CONTINUE 

1203 C COMPUTE HEAT CAPACITY OF GAS MIXTURE AT TGX 

1204 CPGAS=FXN2 (Y, TGX) 

1205 C Limiting Nu and Sh numbers are 2.976. This is taken from page 

1206 C 241 of Convective Heat Transfer by Burmeister, L. C., John 

1207 C Wiley & Sons, (1983) . This assumes a constant temperature 

1208 C profile (Nu) or concentration profile (Sh). 

1209 6669 KC = 2.976 * DAB / CHHD 

1210 H = 2.976 * K / CHHD 

1211 A = AREF * DEXP (-ENGACT/ R * (TREF-TSX) /TSX/TREF) 

1212 THIELE = CMCWT * DSQRT ( A * DENCAT / DABEFF) 

1213 IF ( (THIELE. GT. 100.0) .AND. (ONCE. NE.l) ) THEN 

1214 WRITE (*,*) 'THEILE MODULUS IS > 100. SETTING THIELE=100’ 

1215 WRITE (*,*) 'TO PREVENT NUMERIC OVERFLOW ERROR' 

1216 ONCE=l 

1217 ENDIF 

1218 IF ( (THIELE. GT. 100.0) ) THIELE=100.0 

1219 N = DTANH ( THIELE ) / THIELE 

1220 CONSTl= N*A*CMCWT*DENCAT 

1221 CSX = CGX * ( KC) / (KC+CONST1) 

1222 TSX =TGX + HEATRX*CSX*CONSTl/H 

1223 IF (FIRST. EQ. .TRUE. ) THEN 

1224 A = AREF*DEXP (-ENGACT/ R* (TREF-TSX) /TSX/TREF) 

1225 THIELE = CMCWT * DSQRT ( A * DENCAT / DABEFF) 

1226 N = DTANH ( THIELE ) / THIELE 

1227 CONSTl= N*A*CMCWT*DENCAT 

1228 CSX = CGX * ( KC) / (KC+CONST1) 

1229 TSX =TGX + HEATRX*CSX*CONSTl/H 

1230 IF (SSS.EQ.'I') TSX=TGX 

1231 IF (P.EQ.0) THEN 

1232 XZ=0 

1233 CONVXZ=0 . 0 

1234 CGXZ=CGX 

1235 CSXZ=CSX 

1236 TSXZ=TSX 

1237 ENDIF 

1238 ENDIF 

1239 RE=1000 . 0*CHKD*CHGVZL*DENGAS/VISCTC 

1240 C0N7X = 100.0* (CGXZ-CGX*TGX/TZERO*PZERO/?RESS) /CGXZ 

1241 C HEADER OUTPUT SECTION 

1242 C Print out Reynold's number the first time through the loop. 

1243 C If First = .true, then print out reynold's number, else don't 

1244 C worry about it. A-lso, print out header and results at position X=0 

1245 IF (FIRST. EQ. .TRUE*. ) THEN 

1246 C Output cf selectable parameter to data file 

1247 CALL OUTPUT (1 , MOST, MFDH, MFDV, MCW, NOCHV, NOCHH, MCWT, 

1248 + Y,MVFLZ,TZERO, PZERO, DENCAT, AREF, ENGACT, 

1249 + EMICRO, EMACRO, RMICRO, RMACRO, SSS, SSS1, SSS2, 

1250 + CONV, MONOL, DX, CDI ) 

1251 CALL HEADER ( 0 , RE , CPGAS , DENGAS, CHGVEL, N, DAB, 

1252 + DABEFF, PRESS, Y,DX) 

1253 FIRST=. FALSE. 

1254 WRITE (*,*) 
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1256 

1257 

1258 

1259 

1260 
1261 
1262 

1263 

1264 

1265 

1266 

1267 

1268 

1269 

1270 

1271 

1272 

1273 

1274 

1275 

1276 

1277 

1278 

1279 

1280 
1281 
1282 

1283 

1284 

1285 

1286 

1287 

1288 

1289 

1290 

1291 

1292 

1293 

1294 

1295 

1296 

1297 
1295 
1299 
1330 

1301 

1302 

1303 

1304 

1305 

1306 

1307 

1308 

1309 

1310 

1311 


WRITE (25,*) 

WRITE (*, 9950) 

WRITE (25, 9950) 

WRITE (*,9951) 

WRITE (25, 9951) 

9950 FORMAT (' Distance ', IX, ’ %Conver ' , 2X, ' 02gas ’, 4X, ' 02wall 2X, 

+ ' 02center ' , 3X, ' Tgas ' , 5X, ' Twall ' , 3X, ' DPress ' ) 

9951 FORMAT ( 1 (mm) ’ , 11X, ' (mMOL/L) ' , IX, ' (mMOL/L) ' , IX, ' (mMOL/L) * , 
+1X, * (Kelvin) ' , IX, ' (Kelvin) ' , 3X, ' (kPa) ’ ) 

WRITE (*, * ) 

WRITE (25,*) 

WRITE (*,10) XZ, CONVXZ, CGXZ*lE+09, CSXZ*lE+09, 

+ (CSXZ*1E+09/DCOSH (THIELE) ) , TZERO, TSXZ, (PZERO-PZERO) *101.325 
WRITE (25, 20) XZ, S, CONVXZ, S, CGXZ*lE+09, S, CSXZ*lE+09, S, 

+ (CSXZ*1E+09/DCOSH (THIELE) ), S, TZERO, S, TSXZ, S, 

+ (PZERO-PZERO) *101 . 325 

10 FORMAT (F8 . 3, IX, F8 . 4, IX, F7 . 4, IX, F7 . 4, IX, F7 . 4, 5X, F7 . 3, IX, F8 . 3 
+, IX, E10 . 3) 

20 FORMAT (F8 . 3, Al, F8 . 4, Al, F7 . 4, Al, F7 . 4, Al, F7 . 4, Al, F7 . 3, Al, F8 . 3 
+ , Al, E10 . 3) 

ELSE 

C COMPUTED VALUES OUTPUTTED TO SCREEN AND FILE AT DISPLAY INTERVAL 
C WRITE RESULTS TO SCREEN WHEN REMAINDER OF CONSTA/FACTOR > 0 
C ELSE, JUST CALCULATE VALUES AND DON’T PRINT THEM OUT. 

IF (MOD (CONSTA, FACTOR) .GT. 0.0) GOTO 7777 
IF (SSS1 .EQ.’S') GOTO 7777 
IF (SSS.EQ.’I') TSX=TGX 
WRITE (* , 10 ) X, CONVX, CGX*lE+09, CSX*lE+09, 

+ (CSX*1E+09/DCOSH (THIELE) ) , TGX, TSX, (PRESS-PZERO) *101 . 325 
WRITE (25, 20) X, S, CONVX, S, CGX*lE+09, S, CSX*lE+09, S, 

+ (CSX*1E+09/DCOSH (THIELE) ) , S, TGX, S, TSX, S, 

+ (PRESS-PZERO) *101 . 325 
7777 ENDIF 

C Reset variables each time Runge Kutta Subroutine is called 
7888 X0 (1) = XNEW(l) 

XO (2) = XNEW (2) 

X0 (3) = XNEW (3) 

C Call Runge Kutta Subroutine: 

CALL'RK4 (DX,X,X0,XNEW) 

C PARMETER UPDATE SECTION 
FGX = XNEW (1 ) 

PRESS = XNEW (2) 

IF (SSS.EQ.'I') THEN 
TGX = TZERO 

ELSE 

TGX = XNEW (3) 

ENDIF 

CGX = FGX/CHFLO 
DCGX=CGX-CGXTMP 
X = X + DX 
CGXTMP=CGX 

Y (5) =Y50* (100 . 0-CONVX) /100.0 
Y (4)=Y40-2* (Y50-Y (5) ) 

Y (3) =Y30+2* (Y50-Y (5) ) 

Y (1 ) =Y10 
Y(2)=Y20 
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1312 Y(6)=Y60 

1313 CONST7 = Y (1 ) + Y(2) + Y(3) + Y(4)+ Y(5)+ Y(6) 

1314 DO 33445 JJ=1,6 

1315 Y ( JJ) =Y ( JJ) / CONST7 

1316 33445 CONTINUE 

1317 C ADAPTIVE STEP SIZE SECTION 

1318 IF ( ( (-100 . *DCGX/CGX . GT . . 08) ) .AND . ( JUMP.LT. 1) ) THEN 

1319 DX = DX/5.0 

1320 FACTOR = FACTOR* 5 

1321 IF (DX.LT. IE-4) JUMP=1 

1322 CONSTA=F ACTOR- 1 

1323 WRITE (*,*) 'Reduced step size to DX, ' mm' 

1324 ENDIF 

1325 CONSTA=CONSTA + 1 

1326 C PARAMETER BOUNDARY CHECK 

1327 C Check to make sure that Temperature, Pressure, and Reynold' 

1328 C number are still in valid range. 

1329 C Temperature must be between 200-700 K 

1330 C Pressure must be between 10E-03 and 3.45 bar 


1331 C 

Reynolds number must be 

in laminar region. < 2300 

1332 

IF (TGX.GT. 700.0) 

THEN 

1333 

WRITE (*,*) 


1334 

WRITE (*,*) 1 

Current Gas Temperature = 1 , TGX 

1335 

WRITE (*,*) 


1336 

WRITE (*,*) 1 

Error: 1 

1337 

WRITE (*, * ) 1 

This value exceeds the allowable 1 

1338 

WRITE (*,*) 1 

maximum temperature of 700 K r 

1339 

WRITE (*, *) 1 

Lower the inlet gas temperature . 1 

1340 

WRITE (*,*) 


1341 

WRITE (*, *) 

1 Hit Return When Finished Viewing 1 

1342 

WRITE (*,*) 

1343 

READ (*,621) 

S 

1344 

WRITE (*, * ) 


1345 

WRITE (*,*) 


1346 

GOTO 7564 


1347 

ELSE IF (RE. GT. 2300.0) THEN 

1348 

WRITE (*, *) 


1349 

WRITE (*,*) ' 

Current Reynolds number = 1 , RE 

1350 

WRITE (*, * ) 


1351 

WRITE (*, *) ' 

Error: f 

1352 

WRITE (*,*) ' 

This value exceeds the maximum 1 

1253 

WRITE (*, *) ' 

allowable Reynolds number of 2300. 1 

1354 

WR.ITE (*, *) ' 

Reduce the inlet gas flowrate. 1 

1355 

WRITE (*,*) 

1356 

WRITE (*,*) 

1 Hit Return When Finished Viewing 1 

1357 

WRITE (*,*) 

1258 

READ (*,621) 

S 

1359 

WRITE (*,*) 


1360 

WRITE (*, * ) 


1361 

GOTO 7564 


1362 

ELSE IF ( (TGX.LT. 200 . 0) ) THEN 

1363 

WRITE (*,*) 


1364 

WRITE (*,*)’ 

Current Gas Temperature = ' , TGX 

1365 

WRITE (*,*) 

1366 

WRITE (*,*) 1 

Error: 1 

1367 

WRITE (*, * ) 1 

This value is less than the minimum 

1368 

WRITE (*,*) 1 

allowable Temperature of 200 K. ' 
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1369 

1370 

1371 

1372 

1373 

1374 

1375 

1376 

1377 

1378 

1379 

1380 

1381 

1382 

1383 

1384 

1385 

1386 

1387 

1388 

1389 

1390 

1391 

1392 

1393 

1394 

1395 

1396 

1397 

1398 

1399 

1400 

1401 

1402 

1403 

1404 

1405 

1406 

1407 

1408 

1409 

1410 

1411 

1412 


WRITE (*,*)' Raise the inlet gas temperature.' 

WRITE (*,*) 

WRITE (*,*) 'Hit Return When Finished Viewing' 

WRITE (*,*) 

READ (*, 621) S 
WRITE (*,*) 

WRITE (*,*) 

GOTO 7564 

ELSE IF (PRESS*101. 325. GT. 350.0) THEN 
WRITE (*,*) 

WRITE (*,*) 'Current Pressure* ' ,PRESS*101.325, ' kPa' 
WRITE (*, *) 

WRITE(*,*)' Error:' 

WRITE!*,*) ' This value exceeds the maximum ' 
WRITE!*,*)' allowable Pressure of 350 kPa. ' 
WRITE!*,*)' Reduce the inlet gas pressure.' 

WRITE!*,*) 

WRITE!*,*) 'Hit Return When Finished Viewing' 
WRITE!*,*) 

READ!*, 621) S 
WRITE!*,*) 

WRITE!*,*) 

GOTO 7564 

ELSE IF (PRESS*101. 325.LT. 0.1) THEN 
WRITE!*,*) 

WRITE (*,*) 'Current Pressure* ' ,PRESS*101. 325, ' kPa’ 
WRITE!*,*) 

WRITE!*,*)' Error:' 

WRITE!*,*) ' This value is below the minimum ' 
WRITE!*,*)' allowable Pressure of 0.1 kPa. ' 
WRITE!*,*)' Raise the inlet gas pressure.' 

WRITE!*,*) 

WRITE!*,*) 'Hit Return When Finished Viewing' 
WRITE!*,*) 

READ!*, 621) S 
WRITE!*,*) 

WRITE!*,*) 

GOTO 7564 

ENDIF 

8000 CONTINUE 

q ****** Q r main COMPUTATION LOOP************************************ 


1413 


1414 C OUTPUT OF FINAL PARAMETERS VALUES TO SCREEN AND FILE 

1415 9900 CONTINUE 

1416 16950 CSX = CGX * ( KC) / (KC+CONST1) 

1417 TSX =TGX + HEATRX* CSX* CONST1 / H 

1418 IF (SSS.EQ.'I 1 ) TSX=TGX 

1419 CONVX = 100.0* (CGXZ-CGX*TGX/TZERO*PZERO/PRESS)/CGX2 

1420 WRITE (*,10) X,CONVX,CGX’1E+09,CSX*1E+09, 

1421 + (CSX*1E+09/DCOSH (THIELE) ) ,TGX, TSX, (PRESS-PZERO) *101 . 325 

1422 WRITE (25, 20) X, S, CONVX, S, CGX*lE+09, S, CSX*lE+09, S, 

1423 + ( CSX* 1E+0 9 /DCOSH (THIELE) ) , S, TGX, S, TSX, S, 

1424 + (PRESS-PZERO) *101. 325 

1425 16900 WRITE (25,*) 
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1426 

1427 

1428 

1429 

1430 

1431 

1432 

1433 

1434 

1435 

1436 

1437 

1438 

1439 

1440 

1441 

1442 

1443 

1444 

1445 

1446 

1447 

1448 

1449 

1450 

1451 

1452 

1453 

1454 

1455 

1456 

1457 

1458 

1459 

1460 

1461 

1462 

1463 

1464 

1465 

1466 
146 " 


2^69 

1470 

1471 

1472 

1473 

1474 

1475 

1476 

1477 

1478 

1479 

1480 

1481 

1482 


RE=1000 . 0*CHHD*CHGVEL*DENGAS/VISCTC 
WRITE (*,*) 

CALL HEADER (1 , RE, CPGAS, DENGAS, CHGVEL, N, DAB, 

+ DABEFF, PRESS, Y,DX) 

WRITE ( * , * ) BEEP , BEEP , BEEP 
17000 WRITE (* , * ) 

WRITE (*,*) ' — > Hit Return to Return to Main Menu <==' 

WRITE ( * , * ) 

READ (*, 621) S 
WRITE (*, *) 

WRITE (*,*) 

C RESET OTHER PARAMETERS BACK TO ORIGINAL VALUES. 

7564 MVFLOMVFLZ 

Y (1)=Y10 

Y (2) =Y20 
Y (3) =Y30 

Y (4) =Y40 
Y (5) =Y50 
Y(6)=Y60 
GOTO 9001 

99999 END 

0***************** Qp M&iN ppoGRAM********************************** 


Q************************SUgRQU<pjJJ£ RR4 ******************************* 

C-* Programmer: Seth Daniel Goldblum, UCSD, July 1988 

C-* 

C-* Runge Kutta Subroutine to solve N differential 

C-* Equations simultaneously. 

C-* 

C-* This subroutine relies on a main program to feed it 

C-* the following: 

C-* 

C-* 1) The starting value of the independent variable, P 

C-* 

C-* 2) The step size, DX 

C-* 

C-* The subroutine returns the values of all dependent variables 

C-* at the evaluated point of the independent variable in an array 

C-* Xnew(N) 

C-* . 

C-* The subroutine requires and additional subroutine called func 
C-* which is declared external and contains the functions cf ail 
C-* derivatives to be evaluated simultaneously. 

C-* 


C-* 

C-* 

C-* 

C-* 

C-* 

C-* 

C-* 

C-* 

C-* 

C-* 

C-* 


XO (N) array of N initial values of the functions 

P initial value of the independent variable 

F (N) array of N elements containing the values of 
the functions evaluated at some point. 

DX step size 

XNEW (N) array of N values containing the values of 
the dependent variables 

K1 (N) array of N values containing the values of K1 
used in the Runge - Kutta method 



* 
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1483 

C-* 

K2 (N) array of N values containing the values of K2 


1484 

C-* 

used in the Runge-Kutta method 


1485 

C-* 

K3 (N) array of N values containing the values of K3 


1486 

C-* 

used in the Runge - Kutta method 


1487 

C-* 

K4 (N) array of N values containing the values of K4 


1488 

C-* 

used in the Runge-Kutta method 


1489 

C-* 



1490 


SUBROUTINE RK4 (DX,P,XO,XNEW) 


1491 


REAL* 8 XO (3) , P, F (3) , DX, XNEW(3), K1 (3) , K2(3), K3(3), K4 (3) 


1492 

C-* 

first approx. 


1493 


CALL FUNC (P, X0, F) 


1494 


DO 10 J=l, 3 


1495 


K1 (J) =DX*F(J) 


1496 


XNEW ( J) =X0 ( J) + Kl(J)/2.0 


1497 

10 

CONTINUE 


1498 

C-* 

second approx. 


1499 


CALL FUNC (DX/2 . 0 + P,XNEW,F) 


1500 


DO 20 J=l, 3 


1501 


K2 (J)=F(J)*DX 


1502 


XNEW ( J) =X0 ( J) + K2(J)/2.0 


1503 

20 

CONTINUE 


1504 

c-* 

third approx. 


1505 


CALL FUNC (DX/2. 0 + P, XNEW,F) 


1506 


DO 30 J=l, 3 

! 

1507 


K3(J)=F(J)*DX 

1 

1508 


XNEW ( J) =X0 ( J) + K3 ( J) 


1509 

30 

CONTINUE 


1510 

c-* 

fourth approx. 


1511 


CALL FUNC (DX + P,XNEW,F) 


1512 


DO 40 J=l, 3 


1513 


K4 (J)=F(J) *DX 

i 

1514 


XNEW ( J) =X0 ( J) + (K1(J) + 2*K2 ( J) +2*K3 ( J) +K4(J))/6.0 

i 

1515 

40 

CONTINUE 


1516 


RETURN 


1517 


END 


1518 




1519 




1520 




1521 

£★★★*★ ★★★****************2U£RQU^jNE FUNC********************** *★★★*★★★* 


1522 


SUBROUTINE FUNC(P,X0,F) 


1523 

C-* 

This subroutine is used to store the functions used 


1524 

C-* 

in the runge-xutta method fcr solving them. 

j 

1525 


COMMON /ELK1 / KC , CHSTV, MCH, MOW, CSX, CGX, CHFLO, CPGAS 


1526 


COMMON /SLK2/ VISCTC, CHGVEL, CHHD, K, TSX, DENGAS 

i 

1527 


REAL-8 F (3) ,X0 (3) 


1525 


REAL- 8 KC, CHSTV, MCH, MCW, CSX, CGX, CHFLO, CPGAS 


1529 


PEAL* 8 VISCTC, CHGVEL, CHHD, H, TSX, DENGAS 

1 

1530 


FGX=X0 (1) 


1531 


TGX=X0 (3) 


1532 

C-* 

F (1) is the differential equation describing the change 


1533 

C-* 

in molar flowrate of 02 (moles/s) for a change in length: 


1534 


F (1 ) = -KC-CHSTV-MCH-MCW* ( CGX- CSX) 


1535 

C-* 



1536 

C-* 

F(2) is the differential equation describing the change 

i 

1537 

C-* 

in Pressure for a change in length 

i 

i 

1538 


F (2) = -32 . 0* (VISCTC/1000 . 0) *CHGVEL/CHHD/CHHD/101330 . 0 

i 

1539 

C-* 


i 
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ORIGINAL Paw le 
0F P 0 0 R Quality 


1540 

C-* 

F(3) is the differential equation describing the change 

1541 

C-* 

in temperature for a change in length. 


1542 


F (3) = H* CHSTV*MCH*MCW / CHFLO* (TSX-TGX) / ( DENGAS * CP GAS ) 

1543 

C-* 



1544 


RETURN 


1545 


END 


1546 




1547 




1548 




1549 

lines *************************** 

1550 


Subroutine LINES (K) 


1551 


Integer J,K 


1552 


DO 459 J=l, K 


1553 


WRITE (*,*) 


1554 

459 

CONTINUE 


1555 


Return 


1556 


End 


1557 




1558 




1559 




1560 

£★★★★*★ ★★★★★★★★★★★★★★★★sjjbrqutinE header*************** ************ 

1561 


Subroutine HEADER (FLAG, RE, CPGAS, DENGAS, 

CHGVEL, N, DAB, 

1562 


+ DABEFF, PRESS, Y,DX) 


1563 


Real *8 PRESS, RE, CPGAS , DENGAS , CHGVEL , N, DAB , DABEFF , Y ( 6 ) , DX 

1564 


Integer FLAG 


1565 


If (FLAG.EQ. 0) THEN 


1566 


Write (*,*) ' Initial values' 


1567 


Write (25,*) ' Initial values' 


1568 


ENDIF 


1569 


If (FLAG.EQ.l) THEN 


1570 


Write (*,*) ' Final values' 


1571 


Write (25,*) ' Final values' 


1572 


WRITE (*,*)' Mole fractions: ' 


1573 


WRITE (*, 3233) Y (1) , Y (2) , Y (3) 


1574 


WRITE (*, 3234) Y(4) ,Y(5) ,Y(6) 


1575 


WRITE (25, * ) ' Mole fractions: ' 


1576 


WRITE (25, 3233) Y (1) ,Y (2) , Y (3) 


1577 


WRITE (25, 3234) Y (4) , Y (5) , Y (6) 


1578 

3233 

FORMAT (' He = ' ,F6.4, ' Ar = ' 

,F6 . 4, ' C02 = ' , F6 . 4) 

1579 

3234 

FORMAT ( ' CO = ' ,F6.4, ' 02 = ' 

,F6.4, ' N2 = ' ,F6.4) 

1580 


WRITE (*, 6777) PRESS*101 . 325 


1581 

6777 

FORMAT ( ' Gas Pressure (k?a) : 

' ,FS . 3) 

1582 


WRITE (25, 6777) ?RESS*101 . 325 


1583 


ENDIF 


1584 


WRITS (*,10) RE 


1585 

10 

Fermat ( ' Reynolds number = 

' ,F9. 3) 

1586 


WRITE (*,20) CPGAS 


1587 

20 

FcrmatC Gas Heat Capacity (J/K-g) = 

' ,F9 . 3) 

1588 


WRITE (*,30) DENGAS*1000 


1589 

30 

FormatC Gas Density (g/crr>T3) = 

' , F10 . 7 ) 

1590 


WRITE (*,40) CHGVEL 


1591 

40 

FormatC Gas Velocity (mm/s) = 

' , F9 . 3) 

1592 


WRITE (*,50) N 


1593 

50 

Format ( ' Eff ectivness factor = 

' , F9 . 5) 

1594 


WRITE (*,60) DAB/100.0 


1595 

60 

FormatC Bulk Gas Diffusivity (crre2/s) 

= ' , F9 . 5) 

1596 


WRITE (*,70) DABEFF/100 . 0 
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1597 

1598 

1599 

1600 
1601 
1602 

1603 

1604 

1605 

1606 

1607 

1608 

1609 

1610 
1611 
1612 

1613 

1614 

1615 

1616 

1617 

1618 

1619 

1620 
1621 
1622 

1623 

1624 

1625 

1626 

1627 

1628 

1629 

1630 

1631 

1632 

1633 

1634 

1635 

1636 

1637 

1638 

1640 

1641 

1642 

1643 

1644 

1645 

1646 

1647 

1648 

1649 

1650 

1651 

1652 

1653 


29 


70 

80 


Formate Effective Diffusivity (cm^2/s) = 

WRITE (*,80) DX 

Formate Step size (mm) = 

WRITE (25, 10) RE 

CPGAS 

DENGAS*1000 
CHGVEL 
N 

DAB/100 
DABEFF/100 
DX 


WRITE (25, 20) 
WRITE (25, 30) 
WRITE (25, 40) 
WRITE (25, 50) 
WRITE (25, 60) 
WRITE (25, 70) 
WRITE (25, 80) 
Return 
End 


\F9.6) 
' ,F9. 6) 


0*-*********************2UBROU']’jjq£ TCON ***************************** 
Subroutine TCON(T,TCONT) 

Real *8 T, TCONT (6) 

C Compute thermal cond. (TCONT (i) ) for a given species at temp. T 
C Units on TCON are mW/cM-K. 

C The index of TCONT (i) corresponds to the species. 

C Species : l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2 . 

C The following polynomial equations used to compute the 
C thermal conductivity of each species at temperature T 

C have been taken from page 515 of "Properties of Gases and Liquids" 

C 4th Ed. Reid, Prausnitz, & Sherwood 

TCONT (1) = (3.722E-2 + 3.896E-4*T - 7 . 450E-8*T**2 + 

1 1 . 290E-11*T**3) *10.0 

TCONT (2) = (2.714E-3 + 5.540E-5*T - 2.178E-8*T**2 + 

1 5 . 528E-12*T**3) *10 . 0 

TCONT ( 3 ) = (-7.215E-3 + 8.015E-5*T + 5. 477E-9*T**2 - 
1 1 . 053E-11*T**3) *10 . 0 

TCONT ( 4 ) = (5.067E-4 + 9.125E-5*T - 3. 524E-8*T**2 + 

1 8 . 199E-12*T**3) *10 . 0 

TCONT ( 5 ) = (-3.273E-4 + 9.966E-5*T - 3. 743E-8*T**2 + 

1 S . 732E-12*T**3) *10.0 

TCONT ( 6 ) = (3.919E-4 + 9.816E-5*T - 5. 067E-8*T**2 + 

1 1 . 504E-11*T**3) *10 . 0 

Return 
End 


0-*-***^*r^Tr *-*•*-:* •*★■*■* SUBROUTINE PROP '*'***'* * ★ * ■*-*■«*• * ★ * ■* * •* * -k **■ 

Subroutine PROP (Y, TCONT, VISCT, VISCTC, K) 

Real* 8 AX (6, 6) , VISCTC, K, SUM, MW (6) , Y (6) , TCONT (6) 

Real * 8 CONSTE , CONSTD , VISCT ( 6 ) , LAMBMX, CONST" 

Integer I,J 

C Read in molecular weights MW(I) into array. 

C The index of MW correspond to the species. 

C Species : l=He, 2=Ar, 3=C02 , 4=CO, 5=02, 6=N2 . 

DATA MW (1) ,MW (2) ,MW(3) ,MW(4) ,MW(5) ,MW(6) 

+ / 4.0026,39.948,44.01,28.0105,31.9988,28.0134/ 

C Compute parameter matrix AX(i, j) 

C The indexes of A(i,j) correspond to the species. 

C Species : l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2 . 
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1654 

1655 

1656 

1657 

1658 

1659 

1660 
1661 
1662 

1663 

1664 

1665 

1666 

1667 

1668 

1669 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 

1678 

1679 

1680 
1681 
1682 

1683 

1684 

1685 

1686 

1687 

1688 

1689 

1690 

1691 

1692 

1693 

1694 


| r\ 


1698 

1699 
- 

1 I v j‘J 

1701 

1702 

1703 

1704 

1705 

1706 

1707 

1708 

1709 

1710 


C The AX(i, j) matrix is defined on p. 531 of "Properties of Gases and 
C Liquids", 4TH ed.,Reid, Prausnitz,& Sherwood. Both the bulk gas 
C thermal conductivity and the bulk gas viscosity can be computed 
C using the AX (6, 6) matrix. 

VISCTC =0.0 
LAMBMX=0 . 0 
DO 551 I =1,6 

SUM = 0.0 
DO 651 J = 1, 6 

IF(I.EQ.J) THEN 

AX (I, J) =1.0 

ELSE 

CONSTE = MW (I) /MW ( J) 

CONSTD = VISCT (I) / (CONSTE*VISCT ( J) ) 

CONSTD = DSQRT (CONSTD) 

AX (I, J) = (1+CONSTD*CONSTE**0 . 25) **2/ 

* DSQRT ( 8 . 0* (1 . O+CONSTE ) ) 

END IF 

SUM = SUM + Y (I) *AX (I, J) 

651 CONTINUE 

CONSTF = Y(I) / SUM 

LAMBMX = LAMBMX + CONSTF*TCONT (I) 

VISCTC = VISCTC + CONSTF*VISCT (I ) 

551 CONTINUE 

K = LAMBMX/10000 . 0 

RETURN 

END 


0**********************gU3RQUTiNE VIS ***************************** 
Subroutine VIS (T,VISC, VISCT) 

Real *8 T, VISC (6,2:7) , VISCT (6) ,CONST5 

Integer I 

C Compute viscosity (VISCT (i) ) for a given species at temp. T 
C Units on VISCT are centipoises. 

C The index of VISC corresponds to the species. 

C Species: l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2. 

C CONST5 is a temporary value holding constant 
DO 2223 1=1,6 

CONST5=VISC (I, INT (l+T/100 . 0) ) -VISC (I, INT (T/100 . 0) ) 
CCNST5=CONST5* (T/100 . C-INT (T/100. 0) ) 

VISCT (I) =VISC (I, INT (T/100 . 0) ) + CONST5 
2223 CONTINUE 
RETURN 
END 


£★★*★★**★*★* OUTPUT**' * x * ** * * * * * * * ** * * * * "* * * * * * * * 


Subroutine 

+ 

+ 


Integer 
Real *8 
Real *8 


OUTPUT (FLAG, MOST, MFDH, MFDV, MCW, NOCHV, NOCHH, MCWT, 
Y , MVFLZ , TZERO , PZERO , DENCAT , AREF , ENGACT , 

EMICRO, EMACRO, RMICRO, RMACRO, SSS, SSS1 , SSS2, 

CONV, MONOL, DX, CDI) 

FLAG 

MOST, MFDH, MFDV, MCW, NOCHV, NOCHH, MCWT 
Y (6) , MVFLZ , TZERO, PZERO, DENCAT, AREF , ENGACT 
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1711 

1712 

1713 

1714 

1715 

1716 

1717 

1718 

1719 

1720 

1721 

1722 

1723 

1724 

1725 

1726 

1727 

1728 

1729 

1730 

1731 

1732 

1733 

1734 

1735 

1736 

1737 

1738 

1739 

1740 

1741 

1742 

1743 

1744 

1745 

1746 

1747 

1748 

1749 

1750 

*! 7c;i 


1754 

1755 

1756 

1 7 C'7 

1758 

1759 

1760 

1761 

1762 

1763 

1764 

1765 

1766 

1767 


Real * 8 EMICRO, EMACRO, RMICRO, RMACRO 

Character*! SSS, SSS1, SSS2, S 
Real * 8 CONV, MONOL, DX, CDI 

C Flag=0 => Output parameters to screen for review 
C Flag=l => Output parameters to file LDATA 
IF (FLAG.EQ.l) GOTO 40000 
Write (*, 6000) 

6000 FORMAT ( 2 OX, ’ SELECTABLE PARAMETER SUMMARY • ) 

WRITE (*,*)' Monolith Dimensions (mm) : ' 

WRITE (*, 6001) MCST* 2 , MFDH , MFDV 

6001 FORMAT (' Support wall thickness: ' ,F6.2, 3X, 'Face dimension: 
2F7.2, ’ x ' ,F7.2) 

WRITE (*, 6002) MCW, NOCHV, NOCHH 

6002 FORMAT (' Channel inner dimension: ' ,F6. 2, 3X, 'No. Face channels', 
2’: ' ,F6.2, ' x ',F6.2) 

WRITE (*, 6003) MCWT, 100* (MCW**2/ (MCST*2+MCW) **2) 

6003 FORMAT ( ' Active layer thickness: 1 ,F6. 2, 4x, ' % monol . volume ' , 

+’ open : ' , F6.1) 

WRITE (*,*) 

WRITE (*,*)' Monolith inlet parameters:' 

WRITE (*, 6006) Y(3) , Y (4) , Y (5) ,Y(6) , Y (1) ,Y(2) 

6006 FORMAT (' Gas Composition (mole fraction): C02: ' ,F6.4,2x, 

2 ' CO: \F6.4,2X, '02: ' ,F6.4, /, 38X, 'N2: ' ,F6. 4, 2X, ' He: ',F6.4,2X, 
3'Ar: ’,F6.4) 

WRITE (* , 6007 ) MVFLZ, TZERO, PZERO*101 . 325 

6007 FORMAT ( ' Gas Flowrate (liters/s) :' ,F7. 3, ' Gas Temperature (K) : ', 

2F7.2/ ' Inlet Gas Pressure (kPa) : ',F9.3) 

WRITE (*,*) 

WRITE (*,*)' Catalyst Properties:' 

WRITE (*, 6008) DENCAT 

6008 FORMAT ( ' Catalyst Density (g/mm^3) : ’ ,E10.3) 

WRITE (*,6009) AREF 

6009 FORMAT ( ' Reaction rate constant at 298K(mm A 3 /gcat-s) : ',F8.2) 
WRITE (*, 6010) ENGACT 

6010 FORMAT ( ' Activation energy (J/mol) : ’,F10.2) 

WRITE (*, 6004 ) EMICRO, EMACRO 

6004 FORMAT ( ' Void-fraction as micropores: ' ,F4.2, 

2 ' Void-fraction as macropores: ',F4.2) 

WRITE (*, 6005) RMICRO, RMACRO 

6005 FORMAT (' Avg. micropore radius (nm) : ' ,E7.2, 

2 ' Avg. macronore radius (nm) : ’,E7.2) 

IF (SSS .EQ. ' A' } WRITE (*, 6011) 

6011 FORMAT (' Thermal Operation (adiaoatic/iscthermal) : Adiabatic ' ) 

IF (SSS.EQ. '!' ) WRITE (*, 6012) 

6012 FORMAT (' Thermal Operation (adiabatic/isothermal) : Isothermal') 
WRITE (*,*) 

W r^ITE (*,*)' Computational loop parameters:' 

IF (SSS1.EQ. 'S'j WRITE (*, 6013) 

6013 FORMAT (’ Output file (Full Prcfile/Summary) : Summary') 

IF (SSS1.EQ. 'P' ) WRITE (*, 6014) 

6014 FORMAT ( ' Output file (Full Profile/Summary): Full Profile') 

IF (SSS2.EQ.*'0' ) WRITE (’,6015) CONV 

6015 FORMAT (' Termination on (02 conversion/Length): ', 

2 '02 conversion. %: ', F7.3) 

IF (SSS2.EQ. 'L' ) WRITE (*,6016) MONOL 

6016 FORMAT ( ' Termination on (02 conversion/Length): 

2'Length. Length (mm) : ', F7.3) 


I 


Appendix 4 - program listing 


ORIGINAL PAGE IS 
OF POOR QUALITY 


32 


1768 IF (SSS1.EQ. 'P' ) WRITE (*, 6017) DX, CDI 

1769 6017 FORMAT (' Computation loop step size (mm): ' ,F7.4, 

1770 2' Display every ' , F6.2, ' mm') 

1771 IF (SSS1.EQ. 'S' ) WRITE (*, 6018) DX 

1772 6018 FORMAT (' Computation loop step size (mm): ' ,F7.4) 

1773 WRITE (*,*) 

1774 WRITE (*,6019) 

1775 6019 FORMAT (15X, ' => Hit Return When Finished Viewing <=') 

1776 621 FORMAT (A) 

1777 READ (*,621) S 

1778 WRITE (*,*) 

1779 WRITE (*,*) 

1780 GOTO 80000 

1781 40000 WRITE (*, *) ' *****SEE FILE LDATA FOR RESULTS*****' 

1782 WRITE (25,*) 

1783 WRITE (25, 6000) 

1784 WRITE (25,*)' Monolith Dimensions (mm) : ' 

1785 WRITE (25, 6001) MCST*2,MFDH,MFDV 

1786 WRITE (25, 6002) MCW, NOCHV, NOCHH 

1787 WRITE (25, 6003) MCWT, 100* (MCW**2/ (MCST*2+MCW) **2) 

1788 WRITE (25,*) 

1789 WRITE (25,*)' Monolith inlet parameters:' 

1790 WRITE (25, 6006) Y (3) , Y (4) , Y (5) , Y (6) , Y (1) , Y (2) 

1791 WRITE (25, 6007) MVFLZ, TZERO, PZERO* 101 . 325 

1792 WRITE (25,*) 

1793 WRITE (25,*)' Catalyst Properties:' 

1794 WRITE (25, 6008) DENCAT 

1795 WRITE (25, 6009) AREF 

1796 WRITE (25, 6010) ENGACT 

1797 WRITE (25, 6004 ) EMICRO, EMACRO 

1798 WRITE (25, 6005) RMICRO, RMACRO 


1799 

IF 

(SSS.EQ. 1 

'A') 

WRITE (25, 6011) 


1800 

IF 

(SSS.EQ. 1 

'I') 

WRITE (25, 6012) 


1801 

WRI 

TE (25, *) 




1802 

mi 

TE (25, *) ' 

1 Computational loop parameters 

1803 

J p 

(SSS1.EQ. 

'S') 

WRITE (25, 6013) 


1804 

IF 

(SSS1.EQ. 

' P ' ) 

WRITE (25, 6014) 


1805 

IF 

(SSS2.EQ. 

’O') 

WRITE (25, 6015) 

CONV 

1806 

IF 

(SSS2.EQ. 

’L’) 

WRITE (25, 6016) 

MONOL 

1807 

IF 

(SSS1.EQ. 

ipi ) 

WRITE (25, 6017) 

DX, CDI 

1808 

IF 

(SSS1.EQ. 

'S') 

WRITE (25, 6018) 

DX 

1803 

WRI 

TE (25, *) 




• a* q 

WRI 

TEC*,*) 





1811 80000 RETURN 

1812 END 


1815 

1816 c******^**^****ir-r** FUNCTION FXN1 

1817 FUNCTION FXN1 (TDAE, Y5, PRESS, RMICRO, RMACRO, 

1818 + TSX, EMACRO, EMICRO) 

1819 REAL* 8 TDAB, Y 5, PRESS, RMICRO, RMACRO 

1820 REAL* 8 TSX, EMACRO, EMICRO, DKMICR, DKMACR 

1821 REAL* 8 DMACRC, DMICRO 

1822 C Effective diffusivity is computed using the micro-macro 

1823 C pore model of Wakao and Smith as described on pages 170- 

1824 C 171 of Chemical Reactor Analysis and Design by Froment 
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1825 C and Bishoff 

1826 TDAB= (1-Y5) /TDAB/PKESS 

1827 DKMICR =171482 . 0*RMICRO/1E+7*DSQRT (TSX) 

1828 DKMACR =171482 . 0*RMACRO/1E+7*DSQRT (TSX) 

1829 DMICRO =1 .0/ (1 . 0/TDAB + 1.0/DKMICR) 

1830 DMACRO =1 .0/ (1 . 0/TDAB + 1.0 /DKMACR) 

1831 FXN1 = EMACRO**2*DMACRO + (EMICRO**2*DMICRO* (1.0 + 

1832 * 3. 0*EMACRO) ) / (1.0 - EMACRO) 

1833 RETURN 

1834 END 

1835 

1836 

1837 

1838 c****************** FUNCTION FXN2 ★★★★★★★★**★★★★★★★★★★ 

1839 Function FXN2(Y,TGX) 

1840 C Function computes heat capacity of gas mixture using data 

1841 C from the JANAF tables. 

1842 C MW, CP, Y indexes l=He, 2=Ar, 3=C02, 4=CO, 5=02, 5=02, 6=N2 

1843 Real *8 Y (6) ,TGX,CP (6) ,MW(6) 

1844 C Read in molecular weights MW (I) into array. 

1845 DATA MW (1 ) ,MW (2) ,MW (3) ,MW (4) ,MW (5) ,MW (6) 

1846 + / 4.0026,39.948,44.01,28.0105,31.9988,28.0134/ 

1847 CP (1) =20. 786 

1848 CP (2) =20 . 786 

1849 CP (3) =20 . 712+6. 2501E-02+TGX-1 . 188E-5*TGX*TGX 

1850 + -5. 5773E-08*TGX*TGX*TGX + 4 .1875E-11*TGX*TGX*TGX*TGX 

1851 CP (4) =28. 777 + 5. 5297E-03*TGX -3. 1851E-5*TGX*TGX 

1852 + +7. 094 6E-08*TGX*TGX*TGX -4 . 2501E-11*TGX*TGX*TGX*TGX 

1853 CP (5) =31. 033 - 2. 0261E-02*TGX +6 . 3410E-5*TGX*TGX 

1854 + -4. 9951E-08*TGX*TGX*TGX +9. 1266E-12*TGX*TGX*TGX*TGX 

1855 CP (6) =28 . 567 + 7 . 3978E-03*TGX -3 . 5831E-05*TGX*TGX 

1856 + +6 . 9415E-08*TGX*TGX*TGX -3 . 8499E-11*TGX*TGX*TGX*TGX 

1857 FXN2 = Y(1)*CP(1)/MW(1) + Y (2) *CP (2) /MW (2) + 

1858 + Y (3) *CP (3) /MW (3) + Y (4) *CP (4) /MW (4) + 

1859 + Y (5) *CP (5) /MW (5) + Y (6) *CP (6) /MW(6) 

1860 RETURN 

1861 END 

1862 

1863 

1864 

1865 c**********************FUNCTION FXN3** ************ ******* ★*★★*•*★★★* 

1866 Function FXN3(Y,TGX) 

1867 Real* 8 Y (6) ,D (6) ,TGX, CONST 

1868 C Compute diffusivity of 02 in gas mixture 

1869 C Read in binary diffusion coefficients D(I) into array. 

1870 C The index of D correspond to the species in which 02 is diffusing 

1871 C Note: There is no D(5) defined. These are reference values at 

1872 C 298K, 1 atm, units (mm A 2/s) 

1873 C Species: l=He, 2=Ar, 3=C02, 4=CO, 6=N2. 

1874 DATA D (1 ) , D (2) , D (3) , D (4 ) , D (6) 

1875 + /72. 9,21.2,16.4,15.6,22.5/ 

1876 CONST=0 . 0 

1877 DO 8451 1=1,6 

1878 IF (I.EQ.5) GO TO 8451 

1879 CONST- Y (I) / (D (I) * (TGX/298 . 0) **1 . 5) +CONST 

1880 8451 CONTINUE 

1881 FXN3=CONST 


1 


<r 
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1882 RETURN 

1883 END 

1884 

1885 

1886 

1887 0******************+***fun01’ION FXN4 ********************** ********* 

1888 Function FXN4 (Y) 

1889 Real *8 MW (6) , Y (6) , CONST 

1890 C Compute gas average molecular weight 

1891 C MW, Y indexes l=He, 2=Ar, 3=C02, 4=CO, 5=02, 5=02, 6=N2 

1892 C Read in molecular weights MW (I) into array. 

1893 C The index of MW correspond to the species. 

1894 C Species: l=He, 2=Ar, 3=C02, 4=CO, 5=02, 6=N2. 

1895 DATA MW (1 ) ,MW (2) ,MW (3) ,MW (4) ,MW (5) ,MW(6) 

1896 + / 4.0026,39.948,44.01,28.0105,31.9988,28.0134/ 

1897 CONST=0 . 0 

1898 DO 355 1=1,6 

1899 CONST=MW (I ) *Y (I ) +C0NST 

1900 355 CONTINUE 

1901 FXN4=CONST 

1902 RETURN 

1903 END 

1904 

1905 

1906 c************* END OF LISTING ********************* 

1907 


